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: ABSTRACT 

We investigate higher grading integrable generalizations of the afFine Toda systems, where the flat connec- 
Oh. tions defining the models take values in eigensubspaces of an integral gradation of an affine Kac-Moody 

(-H I algebra, with grades varying from I to —I {I > 1). The corresponding target space possesses nontrivial 

vacua and soliton configurations, which can be interpreted as particles of the theory, on the same footing 



X 



as those associated to fundamental fields. The models can also be formulated by a Hamiltonian reduction 
procedure from the so called two-loop WZNW models. We construct the general solution and show the 
5^ , classes corresponding to the solitons. Some of the particles and solitons become massive when the conformal 

symmetry is spontaneously broken by a mechanism with an intriguing topological character and leading to 
a very simple mass formula. The massive fields associated to non zero grade generators obey field equations 
of the Dirac type and may be regarded as matter fields. A special class of models is remarkable. These 
theories possess a U{1) Noether current which, after a special gauge fixing of the conformal symmetry, is 
proportional to a topological current. This leads to the confinement of the matter field inside the solitons, 
which can be regarded as a one dimensional bag model for QCD. These models are also relevent to the study 
of electron self-localization in (quasi)-one-dimensional electron-phonon systems. 



1 Supported in part by the European Community HCM network # CHRXCT920069. 
^On leave from Instituto de Ffsica Teorica, IFT/UNESP - Sao Paulo - SP - Brasil. 

■^Unite Propre du Centre National de la Recherche Scientifique, associee a I'Ecole Normale Superieure et 
a rUniversite de Paris-Sud. 



1 Introduction 



Integrable theories in low dimensions, besides their intrinsic beauty, have become a very 
important tool in the understanding of basic non-perturbative aspects of physical theories. 
They constitute always a laboratory to test ideas on confinement, quantum physics of solitons 
and many others. In some cases they provide realistic models for very interesting phenomena 
in condensed matter physics, statistical mechanics, and in high energy physics under spe- 
cial kinematical conditions. More recently, they unexpectedly reappeared as describing the 
dependence upon coupling constants of the low-energy effective actions for supersymmetric 
Yang-Mills theories in four dimensions (see e.g. ref.[Q). 

In the present paper we introduce a new class of integrable theories in 1 + 1 dimensions, 
presenting very interesting physical properties, and which we hope, will help understanding 
the role of solitons in quantum field theories. The models generalize the abelian and non 
abelian affine Toda theories, in the sense that they contain matter fields coupled to the 
(gauge) Toda fields. The arising equations are affine specialisations of some general system 
in two dimensions |0 which becomes integrable when associated with a Lie algebra of finite 
growth; see and references therein. They represent affine (non-abelian) extensions of the 
corresponding finite system which, as compared with [Q, contain matter fields. 

We introduce the models through a zero curvature condition, where the flat connections 
take values on a affine Kac-Moody algebra Q endowed with an integral gradation. The 
connection has components not only on the 0, ±1 grades, as for the usual Toda fields, but 
also on eigensubspaces of grades varying from I to — /, with / being a positive integer greater 
than unity. The components of the connection with grades ±1, denoted by E±i, are constant 
(field independent) and play a crucial role in specifying the physical properties of the theory. 
Following 1^, the models are made conformally invariant by the introduction of fields in 
the direction of the central term and grading operator of Q. They can also be obtained by 
Hamiltonian reduction from the so-called two-loop WZNW models [^, ^. 

An initial physical motivation for studying such dynamical systems is the same as the 
one for finite systems |Q. Namely, one describes a nontrivial, not necessarily Riemannian 
target space created by the Toda type fields in the presence of some additional matter fields. 
The latter are related with higher flows of the corresponding fiat connection in the trivial 
holomorphic principal fibre bundle Ai x G i — > Ai, where is a two-dimensional manifold 
and G is an exponential mapping of an affine Lie algebra Q. In the same way as for the finite 
systems, using a relevant specialisation of the Inonii-Wigner contraction it is possible, 
for certain models, to eliminate the back reaction of some matter fields to the Toda type 
fields. As a result, when such a procedure is applied to all matter fields, the latter will simply 
propagate in the field of a given Toda solution. 

Here, in contrast with the finite Lie algebra case, where nonabelian Toda systems lead 
to the exactly solvable conformal systems in the presence of a black hole, the corresponding 
target space possesses nontrivial vacua and sohton configurations, which can be interpreted 
as particles of the theory on the same footing as those associated to the fundamental fields. 

The conformal symmetry is spontaneously broken like in the usual abelian and non 
abelian conformal affine Toda systems generating masses for some particles and solitons 
through a Higgs like mechanism. The masses of the fundamental particles are determined by 
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the eigenvalues of the constant operators E^i, appearing in the flat connection. The masses 
of the solitons have a topological character in the sense that they are determined by the 
asymptotic value of the field in the direction of the central term of ^. It turns out that this 
is also related to the eigenvalues of E±i, and for these reasons we are able to obtain a very 
simple and suggestive mass formula for solitons and particles. 

The massive fields associated to non zero grade generators will be found to satisfy two 
dimensional Dirac type equations. They are interpreted as matter fields, and at first sight 
they are c-number fields, so that they would be dubbed as bosons. However, the issue of 
their statistics can only been solved by considering the corresponding quantum field theory. 
It is well known that in two dimensions the statistics of fields depends upon the coupling 
constant, and perhaps that could be quantized such these matter fields become anticom- 
muting operators. An argument in this direction will indeed be given for a special class of 
models. 

The general solution of the system is constructed following the methods of reference 0, 
based on representation theory of afiine Lie algebras. Some new features appear here, due 
to the higher grade fields, which require more delicate techniques to obtain the expression 
of their general solution. We also use the dressing tranformations |^, |T0|, |Tl|, |12|, as an 
alternative, to construct the solutions in the orbit of the vacuum. The soliton solutions are 
obtained through the so-called solitonic specialization [|1^, Q , see also for the nonabelian 
case. According to that, the one-soliton solutions are determined by choosing the constant 
group element, parametrizing the solutions in the orbit of the vacuum, as an exponentiation 
of an eigenvector of the operators E±i, and the multi-solitons, by taking it as a product of 
such exponentials. 

There is a special class of models which present some remarkable physical properties. 
Any integral gradation of an affine Kac-Moody algebra possesses a period such that the 
eigensubspaces, with grades differing by a multiple of that period, have the same structure. 
For the principal gradation, for instance, that period is equal to the Coxeter number. In 
addition, subspaces of grade n and —n are always isomorphic. By choosing the operators Ei 
and E-i, such that one is the image of the other under such isomorphism, and in addition, 
taking / to be equal to the period associated to the gradation, one obtains models possessing 
a special U{1) Noether current depending only on the matter fields. It is then possible, 
under some circunstances, to choose one solution in each orbit of the conformal group, such 
that for these solutions, that U{1) current is equal to a topological current depending only 
on the (gauge) zero grade fields. The submodel obtained by such special gauge fixing of 
the conformal symmetry, presents some very interesting properties due to this equivalence. 
We show for instance, in the case of a model associated to s/(2), that the matter fields get 
confined inside the solitons. In addition, the masses of solitons and particles are shown 
to be proportional to their U{1) charges, in a manner very similar to what occurs in four 
dimensional gauge theories with Higgs in the adjoint representation and in the BPS limit. We 
believe that this equality between topolgical and Noether currents will play an important 
role in the understanding of the quantum theory of the solitons. We also point out that 
such type of models are related to several interesting phenomena in (quasi)-one-dimensional 
electron-phonon physical systems. 

The paper is organized as follows. In Section 2 we define the models, presenting their 
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equations of motion and discussing their symmetries. The general solution, in terms of 
highest weight representations is worked out in detail in Section 3 and illustrated with 
many examples. The holomorphic factorizable representation of the general solution given 



in Section 3, admits a remarkable specialization [|13|, which is used in Section 4 for the 
calculation of the soliton solutions. In Section 5, we perform the related dressing procedure, 
which gives all solutions in the orbit of the vacuum, including the solitonic ones. The 
dressing method gives further insight into the soliton properties and it is very useful for the 
applications which follow. In Section 6, an alternative formulation of our system as two-loop 
WZNW model, provides an improved energy-momentum tensor obtained via the Sugawara 
construction [0, |T^. The masses of solitons and particles are calculated in Section 7 with 
the help of that tensor, through the spontaneous breakdown of the conformal symmetry by a 
Higgs like mechanism. In Section 8 we discuss the physical properties of the higher grading 
fields establishing that the massive ones satisfy Dirac like equations and discussing some 
of their peculiarities, related to parity and complex conjugation. Section 9 is devoted to a 
specially interesting class of models possessing the U{1) Noether current mentioned above. 
In Section 10, we treat in great detail two models associated to the principal gradation of 
sl{2Y^\ discussing their physical applications. Section 11 is devoted to the conclusions and 
perspectives for future investigations. 



2 Formulation of the System 

Consider an untwisted affine Kac-Moody algebra Q endowed with an integral gradation 
Q = 0„eZ denote 

^+ = 0a„, g- = ^g^n. (2.1) 

7i>0 n>0 

Notice that by an affine Lie algebra we mean a loop algebra corresponding to a finite 
dimensional simple Lie algebra Q of rank r, extended by the center C and the derivation 
D. According to ||15|, integral gradations of Q are labelled by a set of co-prime integers 
s = (so. Si, . . . Sr), and the grading operators are given by 

Q, = H, + N,D~^Tt {H,f C . (2.2) 

Here 

H,^j2saX:-H\ N^^j^s^mf, ^ = ^ m^a, , mt = l; (2.3) 

a=l 1=0 a=l 

is an element of the Cartan subalgebra of Q; aa, a = 1, 2, ... r, are its simple roots; i/j is 
its maximal root; the integers in expansion ip = J2a=i ^t'^a] and are the fundamental 
co-weights satisfying the relation ■ = 6ab- 

Let be a two dimensional manifold with local coordinates x+ and a;_; ^ be an affine Lie 
algebra corresponding to a finite dimensional complex simple Lie algebra Q with the Lie group 
G; ^ be a fiat connection in the trivial holomorphic principal fibre bundle Ai x G i — > A4. 
Specify the connection in such a way that its (1, 0)-component takes values in the subspaces 
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®n=o^+n, and (0, l)-component takes values in 0|^=o^-"; with / being a fixed positive 
integer. In other words, up to a relevant gauge tranformation, these components, satisfying 
the zero curvature condition 

d+A^ - d-A+ + [A+ , A_] = 0, (2.4) 

are of the form 

A+ = -BF^ B-\ A^ = -d-B + (2.5) 

Here i? is a mapping from A4 to the Lie group Go with the Lie algebra Qo; are mappings 
to 0^=1 Q±n of the form 

F+ = i?,+ ^F+, F- = i?_,+ ^F-, (2.6) 

m=l m=l 

with E±i being some fixed elements of Q±i] and F^, 1 < m < / — 1, take values in Q±m- 
Substituting the gauge potentials ( p.5|) into (|2.4| ), one gets the equations of motion 



i-i 



d+{d^BB-') = [E^i , B El B-'] + Y.[F- , B F^ B~'] , (2.7) 

n=l 

= [Ei,B-'Fr.^B]+ Y: [F;1+^,B-'F-B], (2.8) 

n=l 

d^F- = -[E^i,BF,t,r.B^']- Y: [F-^m,BF^B-']. (2.9) 



n=l 



Note that a consideration of the systems generated by the flat connection with the compo- 
nents A± taking values in the subspaces ©I^Lo^in where l± can be different positive integers, 
follows completely the same line. The systems of that type, even without appealing to their 
explicit formulation, are integrable, whenever the corresponding Lie algebra is finite dimen- 
sional or affine; see p|, and also [0, [Q. We restrict ourselves here only to the case with 

= = I] however, other (asymmetric) possibilities can be attractive as well, and their 
investigation follows the same arguments as presented here 0. 

Since Qg and C are in Qq, we parametrise B as 

B = be"'^^e''^ , (2.10) 

where 6 is a mapping to Gq, the subgroup of Gq generated by all elements of Qo other than 
Qs and C. The fields rj and u correspond to the extension of the loop algebra, and, as we 



"^Finishing the present paper, we were acquainted with paper 1 19 where the authors have studied some 
special non-left -right symmetric, as they called heterotic conformal Toda system, corresponding to the case 
of the series Ar endowed with the principal gradation, and a choice when = 2, Z_ = 1. (In fact, their 
consideration is pretty valid for an arbitrary finite dimensional Lie algebra, and, with a minor modification, 
for an affine Lie algebra.) They also used some additional restrictions, in particular that, in our notations, 
E2 = EjElE°a+i)l Pi = [^2, Ff], E,2 = 0, Ff = E°_^, where Ff is a mapping to G^,. 
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will show below, are responsible for making the system conformally invariant |^ 0]. Clearly, 
the order of the three factors in (|2.10|) is irrelevant, since they commute. In addition, we 



will use a special basis for the generators of Go such that they are all orthogonal to Qs and 
C. From (|2.2| ) one observes that the generators of Qq are, besides C and Qs, the elements 
H^, a = 1,2, .. .r, of the Cartan subalgebra, and step operators E^^ and such that 
J2a=i Sa^a ' ct = 0, and J2a=i ^a^a ' P = ^s- There can be no step operators E!^, with | n |> 1, 
as explained in appendix C of ref. 0. Therefore, shifting the Cartan elements as 

^° = - ^ Tr {h^ H^) C=H'^- C, (2.11) 



one gets 



Tr {C) = Tr (c H^) = Tr {qI) = Tr (Q, H^) = 0, Tr (Q, C) = iV„ 

Tr H^) = Tr = 4a, ■ a,/alal ^ r/,,, (2.12) 

for all a, 6 = 1 . . . , r. Here we have used = 2aa ■ H^/ a^, Tr (x ■ y ■ H^) = x ■ y, and 
Tr (C Z^) = 1. For more detail of such a special basis, see appendix C of ref. p. 



Substituting (|2.1CI|) into the equations of motion (|2.7|)-(|2T9|), one has 



i-i 

d+(d-bb-') + d+d-iyC = e'''[E^i,bEib-^] + J2e''^[F~,bF+b-'], (2.13) 



n=l 

l—m—l 



= e^'--^^ [E, , b-' Fr_^b] + Y: e-^i^^+n , , (2.14) 



n=l 

l—m—l 



d^F- = -e(^-m)v^E_,,bF,t,^b-']- E e-^[F~^^,bF;lb-'], (2.15) 



n=l 



d+d^riQ, = 0, (2.16) 

where the last equation is a consenquence of the fact that D, and hence Qs, can not be 
obtained as the Lie bracket of any two elements of Q. 

Let us discuss briefly the meaning of system (2.13)-(2.16). First, it is clear that for I = 1, 
when all the mappings F^, n > 1, are absent in the game, these equations coincide with 
the standard conformally affine Toda system. Second, let us conventionally call the fields, 
parametrising the mapping b, the Toda type fields; and the fields, entering a parametrisation 
of the mappings F^, the matter fields coupled, in accordance with equations (2.13)-(2.16), to 
the Toda type fields; the reason for that will be clear from the following observation. Namely, 
using a relevant specialisation of the Inonii-Wigner contraction [§], one can bring to zero 
the back reaction to the Toda type fields for some or all of the matter fields. Recall that 
under the Inonii-Wigner contraction of a simple Lie algebra Q, its elements are multiplied by 
constant parameters, some of which tend to zero in a consistent way. Roughly speaking, with 
this procedure one multiplies by a contraction parameter, say k, the elements of a subspace 
V C Q, complementary to some subalgebra H of Q; V ^ kV = V^'^\ H ^ H. Then it is 
clear that there exists a limit k ^ 0, when the corresponding algebra Q^'^^ contains P^^-* as an 
ideal; in other words, Q becomes the semi-direct sum Q^^^ of Ti. and V^^\ and, hence, we end 
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up with a non-semisimple Lie algebra. In particular, if a simple Lie algebra Q is endowed 
with a Z-gradation, supply the grading subspaces Qn, corresponding to the subspace V, with 
the parameter k, and then tend it to zero in a consistent way, having in mind the grading 
property [Qm, Gn] C Gm+n- The same scheme takes place for the affine algebras. Then, we 
can eliminate contributions coming from some appropriate mappings in (2.13)-(2.15), 
and, of course, have a possibility to arrive at the case when equation (2.13) does not contain 
the sum of the last / — 1 terms in the r.h.s. As a result, we come to an equation which looks 
similar to those for the standard affine Toda system, however, in general, with a different 
meaning of the elements E±i which belong here to the subspaces Q±i. Evidently, there are 
many other meaningful possibilities. Note that an analogous Inonii-Wigner procedure, albeit 
for the corresponding representations of the algebras, can be applied to obtain the solution of 
the contracted systems, starting from the general solution to system (2.13)-(2.16) obtained 
in the next section. 

The structure of the vacuum of the system ( p.l3|) -( p.l6|) is rather complicated. We will 
discuss some aspects of it below. However, there is a simple condition that guarantees the 
existence of static (vacuum) solutions. If the elements E±i satisfy the relation 

[El , E_i] = /3C , where /3 = ^ Tr {Ei E.i), (2.17) 

then 

6=1, F^ = 0, r? = 0, p = -l3x+x-, (2.18) 
is a (vacuum) solution of (|2.13|) -( p.l6|) . 



Another possibility for vacuum solutions arises when E±i, / > 1, belong to a Heisenberg 
subalgebra of Q, see [jl5|, |20| , 

[Em. En] = Tr (E^^-a/) Mdu+NfiC, (2.19) 

where M,N belong to some (infinite) subset of the integer numbers 2. In such cases 
one has that 

6 = 1, r] = 0, F^, = ci,E±M, F^ = 0, if m ^ , iy = -nx+x_, (2.20) 
is a solution of ( ^.13D -( p.l6D with being constants, and 



i-i 

n = P+J2TT {EmE_m) Met, c^. (2.21) 

A/=l 



Obviously, the system ( 2.13| )- (|2.16 ) may have many more vacuum solutions besides ( p. 18 



and ( p.20| ). However, the condition ( p.l7| ) guarantees the existence of at least one vacuum 



solution. Such a fact, as we will see below, favors the existence of soliton solutions. 

The models introduced above are completely characterised by the data {Q,Qs,l,E±i}] 
and we have a quite large class of systems with physical properties crucially depending on a 
choice of those data. 

Equations ( |2.13| ) - ( p.l6|) are invariant under the conformal transformation 



x+^fix+), x^^gix.), (2.22) 
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with / and g being analytic functions; and with the fields transforming as 

b{x^ , x^) = b{x+ , x^) , (2.23) 

^-H^+,x.) _^ g-i>(x+,S_) _ ^2.24) 
g-r,(x+,x_) ^ g-rj(i+,x_) ^ (^^y// (^^yAg-r,(x+,x_) ^ (^2.25) 

^ F+(i;+,5„) = (/')-'+™/'F+(x+,x_), (2.26) 
F-(x+ , X-) F-{£+ , X-) = {gT^^""^^ F-{x+ , x.) , (2.27) 

where the conformal weight 6, associated to e'", is arbitrary. 

Notice that the Lorentz transformation x± X^^x± is obtained from (|2.22|) by taking 
f{x+) = X+/X and g{x-) = Ax_. Therefore, from (|2.23|) - (|2.27|) we get the fields transforming 
as 

Equations ( |2.13| )-( ^.16| ) are also invariant under the transformations 

6(a;+,x_) — s> /iL (x_) 6 (x+ , x_.) /li? (x+) , (2.29) 
F+(x+,x_) ^ /i^i(x+)i^+(x+,x_)/iK(x+), (2.30) 
F-(a;+,x_) ^ /iL (a;_) F" (x+ , x_) /i^^ (x_) , (2.31) 

where hi (a;_) and (x+) are elements of subgroups TIq and of Go, respectively, satis- 
fying the conditions 

hR{x+) Eih~^\x+)=Ei, h-^^x^) E^ihL{x_) = E_i. (2.32) 

The left and right gauge transformations commute, and so the gauge group is TYq ® TCq. 
Whenever TYq and Hq have a set of common generators, we get an important subgroup of 
the gauge group, namely TCd = 'Hq fl TIq. These are global gauge transformations, where 
the fields are transformed under conjugation {h^ = = hr, = const.), 

h ho hh-^^ , F^^hD F^h~^\ (2.33) 

and E±i = hn E±ih]^^. We discuss the relevance of these transformations below. 



3 General Solution 

The fiat connection ( |2.5| ) and the equations of motion themselves are written in a completely 
analogous way as those for the corresponding finite dimensional systems and, as we have 
already mentioned in the Introduction, they represent a specialisation of a more general 
scheme discussed in on the basis of some original papers refered therein. To obtain the 
general solution, we have to look at the problem of a similarity between the structure of the 
highest weight representation spaces of Q and Q, annihilated by the action of the subspaces 
Q+m and Q+m, respectively. To clarify this point, let us recall, following some results 
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concerning the corresponding finite dimensional case, however, on some more general level, 
applicable for the affine systems. 

First, give some definitions and notations. The flat connection A in question, is repre- 
sented in the gradient form, 

A± = g-'d^g, (3.1) 
with g : Ai i — > G; and we write A± with the help of the modified Gauss decomposition for 

g, 

g = fJ'-i^+lo- and g = fi+i^.-fo+, (3.2) 

respectively, in accordance with the Lie algebra decomposition Q = Q- Q) Go Q) G+; 7o± '■ 
M. I — > Go, lJi'±,v± '■ M. I — > G±. The grading conditions provide the holomorphic property 
of /i-t, namely that they satisfy the initial value problem 

d±fi±{z±) = ^i±{z±)£±{z±), (3.3) 

where 

M 

£±{z^) = ^ti^^y^ ^^($^) = E (3-4) 

™=1 aeA+ 

with arbitrary functions ^^"^{z±) determining the general solution to the system under 
consideration; A+ is the set of the positive roots of ^ = Y.m&z Qm corresponding to the root 
vectors Xa in the subspace Qm- Denote by \h >*^'"^ and the dual, ^'"^ < h'\, to \h >^^\ 
the basis vectors in a representation space, annihilated by all the subspaces Qj^n and 
n > m, respectively. Then the general solution for the Toda type fields contained in B, is 
expressed in terms of 

< h'\{j+)-'fi-'fi^%\h >«= « < h'\B-'\h >(i); (3.5) 

cf. while the others, say matter fields from V^, are determined via the matrix elements 

« < h'\go'fi-'fi^\h >(™), and ('"^ < /iV;V-^o^> >^'^ • (3.6) 

Here = 7o+7o-^ 70^ = 70^(^=1:) • ' — ^ '^o are arbitrary mappings; and B~^ = 
{'-fQ {x+))^^go'yo (x-). The whole set of arbitrary functions which determine the general 
solution to the problem, consists of those parametrising the mappings 7^, and the functions 
a e A+, 1 < m < M - 1; while are expressed in terms of corresponding parame- 
ters of 7^. The mappings are determined via the matrix elements ( p.6|) in a similar way 
to the construction discussed very briefly in [Q. Here we give this construction with more 
details. 

Namely, comparing decompositions ( |2.5| ) with those obtained by substitution of ( p.2| ) in 
( p.l| ) with account of the holomorphic property of fi±, one sees that 

= T{l^)-'B^\v^'d^v^)B^'^l (3.7) 



Consider, for example, the first matrix element in p.6| ), which can be identically rewritten 
as 

< /i'|^o"V;V-|/^ >^™^= < h'\{%\Sgo)u^^\h « < h'\u-^\h . (3.8) 
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Then, differentiating ( ^.81 ) over with the help of ( |3.7| ), we have 

< /i'|^o"V;V-|/i >^"^] = < h'\{BF+B-^)u-^\h (3.9) 
where B = B'-fQ. Now, with m = 2, relation (|3.9|) takes the form 

< h'\g^'fi-'fi^\h >(2)] = (1) < h'\BF+B~'\h >(2), (3.10) 

which determines the mapping F^ . To make the next step and calculate F2 , recall that the 
mapping entering equation (|3^ ), is expressed in terms of £ = BF^B~^ as the Volterra 
type decomposition 

CO „ 

^;' = E / , dyC{y')C{y')---C{yn, (3.11) 

where n"(x+) is the simplex in R" given by Q'"'{x^) = {y G R" '■ a < y^ < y^~^ < ■ ■ • < 
y^ < a is some constant determining the initial value problem; z/+(a) is taken to be 
unity, with an appropriate choice of the initial value problem for n±. (Note that here, of 
course, the mappings C{y) depend also on the variable X-.) Then, with account of formula 
( p.ll| ), it follows from (|3.9| ) that for m = 3 there is a relation 

W < h'\BF+B~^ p dy{BF+B-^)\h >(=^) < h'\BF+B'^\h >(3), 

and hence, knowing already F^ , we obtain an expression for F^. Continuing the procedure, 
we determine all the mappings by the recurrent formula 

< h'\BF+B-^\h >(™+i)= < /i'|^o~V;V- - (3.13) 
f 5: / dyC,,{y')C,,{y')---C,Ayn\h>^-^% 

where Ck^ = BF^^B^^ with integers ki>l entering the sum in (|3.13|) only for Yl,l=i ki = m 
and only once if ki = ki^i. An analogous scheme allows to calculate the mappings F~. 
Note in this context that, in comparison with the method for constructing the Toda type 
fields containing in S, which is, of course, a direct repeatition of those in the proposed 
construction of the matter fields F^ represents a novel feature of the integration scheme. 



The matrix elements ( p.5|) and (|3.6| ) realise, in general, a matrix version of the tau- 



function for the standard Toda system. As it was already noted in from the point of 
view of the Lie algebra representation theory, they are closely related to the Shapovalov form 
defined on a Lie algebra Q. Such a relationship seems to be quite general and natural; it 
takes place for a wide class of nonlinear integrable systems, including, in particular, abelian 
and nonabelian (conformal and affine) Toda systems associated with the simple Lie algebras 



Q. First recall some known definitions, see e.g. |]22[, concerning the Shapovalov form. 

Let p be a Cartan subalgebra of p* be an algebra dual to p; U{Q) be the universal 
enveloping algebra for Q. The Shapovalov form defines the linear mapping U {Q) ®cU{Q) 1 — ^ 
U{p) and is realised as a bilinear form {x^y)^ for any two elements x,y & U{Q). Here 
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X — > is the Chevalley involution for = x', and the hermitean Chevalley involution 
for x^ = X*; the subscript "0" means the projection of U{Q) on U{p) which is parallel to 
Q-U{Q) + U{Q)Q^. Note that the given definition is naturally extended for the case of the 
algebra U\Q) = U{Q) ^u{p) R{p*), where R{p*) is the algebra of rational functions over 
p*. It is very important that the form (x^?/)o is degenerated on the left ideal U'{Q)Q+, and 
is not degenerated on the subalgebra U{Q_)] and hence is not degenerated on the space 
U' {G) /U' {G)G+ which is a rational span of the corresponding Verma module. 

One can get convinced that the matrix elements (p.5| ) and ( p.6[ ) are nothing but the 
Shapovalov type forms (x^y)o for some special two elements x,y E Uh{G-) of the Lie algebra 
G-, with the coefficients being holomorphic functions. And, moreover, since this consideration 
is valid for the affine case as well, the general solution to equations ( p.7|) - ( p.9|) which are the 
conformally affine analogues of the corresponding conformal system [Q, also can be written 
in terms of these forms. The reason is that, in accordance with the general construction, see 
and references therein, the grading conditions, realised in the form of decomposition (|2.5| ) 
for the connection A, provide the holomorphic factorisability of the solution; a differential 
geometry formulation of this fact is given in [ p3[] . 

An explicit formulation of solutions like ( p.5| ) and ( |3.6| ) as the series (infinite for an affine 
case, while absolutely convergent in accordance with the arguments given in and references 
therein) in arbitrary functions ^"^"^ can be done, when it is needed, by the following purely 
algebraic scheme, standard for the Lie algebra representations theory. Construct a basis in 
U{G) with the help of the monomials = Xm„ ■ ■ ■ in the basis elements Xm = Xmi^m) 
of G- In particular, let X-tm, X_m = X+m, be such a basis in U{G±) with the weight /im- 
Then the elements Xj^^X^^ generate a basis of U'{G) over -R(p*), and this procedure gives 
for any weight /i G p* a vector space F^{G) of all formal series Z]m,n Cm,nX+mX_n with 
Cm,n G where the sum runs over all monomials of the weight /i = /im — A^n- The 

subspaces F^{G) are in turn the subspaces of the algebra F{G) graded by the weights /i. 
Finally, one needs to transform the elements X+mX_n entering (|3.5|) and ( |3.6|) to the series 
of the monomials X_XoX+, with Xq e U{p), by using the commutation relations of the 
algebra G- Note that for the abelian systems corresponding to the principal gradation of G-, 
the most suitable choice is the Verma type basis when the monomials Xm are constructed in 
terms of the Chevalley generators satisfying the defining relations. For other gradations, an 
adequate basis has not been discovered yet, and this is why the known explicit expressions 
for the general solutions of the nonabelian Toda type systems are written in a rather 
complicated form, though their holomorphic factorisability in terms of the form like (|3.5|) , 
related in this case with the generalised Verma modules |Q, is quite clear. 

Let us now show how the structure of the matrix elements (|3.5|) and ( p.6|) is made concrete 
for the case of the principal gradation. Here it is suitable to take as the basis vectors \h >^^^ 
the highest weight vectors |i>, l<z<r, of the fundamental representations of ^. It seems 
more illustrative to discuss the formulation in question for the finite Toda systems; for the 
corresponding affine deformations, due to the similarity between the structure of the highest 
weight representation spaces of G and G-, it is practically the same. 

Since for the principal gradation of G-, the subalgebra Go is abelian, parametrise the 
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mapping and 7o'(a;±) as 

where (0, //O) = ELi 0*^°- 

To write down expressions ( p.5|) and ( p^) in an explicit way, we use the Verma type basis 
for the spaces of the fundamental representations of Q. There, the highest vector \i > of the 
i-th fundamental representation satisfies the relations 

E^\i >= 0; H^\t >= 6^J\^ >; E^_j\t >= 0, j; (3.14) 

and the basis we are looking for, are constructed only with the help of the Chevalley gener- 
ators E'^j] namely we choose as the basis vectors those of the set 

\jm...3i;t> = ^°,.^---£;°,.j2>, l<J<r, 0<m<iV,-l; (3.15) 
\jo...ji;i> = \i>, 

with nonzero norm; Ni is the dimension of the representation. There take place the following 
useful formulas: 

E°j \im ■■■ji;i> = \jjm ...ji;i>, 

m 

Hj\jm- ■ -Jui > = i^ij ■ -j^'^^ >^ (3-16) 

9=1 

m q—l 

Elj\jm...ji;i> = ^iU - )\jm---jq---ji;i>, 

q=l n=l 

which are evident in virtue of the defining relations for the Cartan and Chevalley elements, 
namely, 

[H^,H^]=0, [HlE%] = ±k,.E%, [E%E%] = 5.,H^, (3.17) 

with k being the Cartan matrix of Q. Here jq means the absence of the root vector E^j^ 
in the corresponding formula for the basis vector \jm ■ ■ >■ If the norm of the vector 
\jjm ■ ■ ■ > or \jm . . . jq . . . > is equal to zero, the corresponding term in the r.h.s. of 
formulas ( 3.16 ) is naturally absent. 

For example, let us illustrate the structure of the fundamental representation space by 
an example of the simple Lie algebras of rank 2, i.e., the algebras A2 = s/(3, C), B2 = 
0(5, C), C2 = sp{4:, C), and G2- Namely, the basis vectors of the 1-st and 2-nd fundamental 
representations of these algebras respectively are 

A2: 







|2i;i) 


|2), 


|2;2), 


|12;2) 



Bo: 



|1), |1;1), |21;1), |221;1), |1221;1); 
|2), |2;2), |12;2), |212;2); 
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|2i;i), 


|121;1); 


|2), 


|2;2), 


|12;2), 


|112;2), |2112;2) 



G2: 

|1), |1;1), |21;1), |121;1), |1121;1), |21121; 1), |121121;1); 
|2), |2;2), |12;2), |112;2), |1112;2), |2112;2), |21112;2); 
1221112; 2), 1121112; 2), |2121112;2), |1221112;2); 
112121112; 2), 1112121112; 2), |2112121112; 2). 

In the case under consideration here, one can take \jm-i ■ ■ -ji'-ii > as the basis vectors 
1/1 >(™\ Then, using the relation Hj\i >= 6.ij\i >, expression (|3.5|) is written in the form 

g-<^. = e^^-^-" < i|/i;V-K >; (3.18) 

while the first formula in ( |3.6| ) reads as 

W < /i'|7o-6(7o+)"V;V-|/i >^'"^= e^^^'^^-'^^ < ^|^;V-|jm-i ■■■Jl■,^>■ 
Finally, the recurrent formula ( p.l3| ) is reduced to the expression 

<^\F;^\J„,...J^;^>= (3.19) 
eJ27=iT.lik,^i(^--^)id+{e^^+^^-'f'^ < i|/i;V-|jm . . . Ji; ^ > - 

m „ 

/ < ^1^1(2/^)^2(1/^) ••■ A. (y") I jm- • >}• 

Here, with the decomposition = J2aeA+, fa^+a^ mappings read as 

where the product in the exponential is defined by the commutation relation [(0, H^), -E+q,] = 
(a, (l))E^^. All the matrix elements of type 

< ^l/W^V-lim • • • > and < i\Cki{y^)Ck2{y'^) ■ ■ ■ CkXy')\3m ■ • • > 

can be calculated in the explicit form of finite polynomials, using the defining relations ( p.lT] ) 
and the properties (|3.14|) - (|3.15|) of the basis vectors. Technically, it is the same task as for 
the standard Toda systems, see 0. 

In particular, for m = 1, since = J2j fj^E'jy has 



< ^\F,+ \^■^>= ft = ei:Li^-'(^'-^)'g+{ ^ ^ }. (3.20) 



< i\ 


/i+V-l 


i;i> 




\i > 
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For m = 2 one can represent the mapping F2 as 

where the matrix Aj^j^ is obtained from 25j^j^ — kjj^ by replacing all nonzero entries by 

/+j2(x+,a;_) = -/+j-^(a;+,x_), so that 

< i\F^\j2jui >= 6ij^{26ij2 - kij^)Aij^. 
Remembering that Ci = J2j e^^''^"'^ fj'Ej, formula (|3.19| ) for m = 2 is read as 

^ + <^|/i;V-K> 

+ h.e^'^'i'-^'^^ft p dye^'^'i'-^'^^^y^f^iy)} (3.21) 

+ hj 77^1 p 0+ -^^rr r. }, « t J- 

Note that the matrix element < i\fi^^fi_\ji; i > also vanishes when kij = 0. 

To obtain explicit solution for the coefficient functions entering the mappings for 
higher values of m, rewrite the above given decomposition as 

m 

K= E /^...™^f::_,...' E e A+, (3.22) 

n---«m S=l 

with tti being the simple roots of Q; thereof 

A,(y^)= E e ^u^^nEU ^ , (3.23) 

1 ■ *i 



A; ■ -1- 

where I]sLi ^ ^t,- Introduce now the notation 



with which 

<^l^mlim---Ji;«>= E fi <i\Ela\3m---3ui> = 



A + 
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and 



< i\Ck^{y^)Ck^{y'^) ■ ■ ■ CkXy')\jra ■ ■ ■ > 



s ■^^^ 



Here the sum over /^^ means the sum over all indices Hence, the recurrent 

formula ( p.l9| ) takes the form 

ft-i = (3-25) 



<ii|/i+V-|jm---ji;« > 



Here the matrix elements ( p.24| ) can be calculated explicitly. For this goal one needs to 
rewrite the root vectors E'^ as monomials over the Chevalley generators, and then use the cor- 
responding formula from for the matrix elements of type < ■ ■ ■ E^^^E^j^ ■ ■ ■ E^j^ \i > 
written in terms of the Cartan matrix. 



4 Solitons from the general solution 



To obtain soliton solutions to the system (|2.7|) -(2.9) from the general solutions described in 
the previous section, we use the method suggested in for the abelian affine (periodic) 

and 



Toda system, and then discussed in [14 , |25 



for their nonabelian generalisations. 



In the cases when E±i, I > 1, live in a Heisenberg subalgebra of Q, see ( p.l9|) , consider the 
initial value problem (p.3|) with S± in (p.4|) taken to be £±, which are 



E. 



±1 



i-i 

E' 

N=l 



'N 



E±N , and so [S+ , S-] = QC; 



(4.1) 



where c% and f2 were introduced in ( |2.20| ) and ( |2.21|) , respectively. In other words, choose 
all arbitrary functions ^"^"^ there to be zero except those standing in the direction of the 
Heisenberg subalgebra generators E±m and E±i. If E±i do not lie in a Heisenberg subalgebra. 



(4.2) 



take £± in ( |3.4|) to be just E±i (c„ = 0). Then the mappings fi± are 



where /i° are some fixed mappings independent on the local coordinates x±. It is quite clear 
that then one ends up with the following expressions for the matrix elements ( |3.5[ 
determining, with such a choice of the mappings, a particular solution to system (|2.7|)-(|2i 



,(/) = M < h'\ 



\h 



(4.3) 
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with /i° = being a fixed mapping J\4 i — > G. Next step in the calculation is to 

remove the dependence on the elements S± in ( |4.3| ), which can be easily done by replacing 
the exponential with S+ {S-) to the extreme right (left) position where it gives unity under 
the action on the corresponding basis state which is annihilated by {S-). Note that this 
procedure is, of course, as it was mentioned in the previous section, common for all methods 
for obtaining an explicit series for the matrix elements under consideration, and is traditional 



in the representation theory, see e.g. [^. Now we make use of the second ingredient of 
the scheme which consists in an assumption that the mapping /i" is an eigenvector with 



respect to £+ and Finally, following [|13[, we write product of the exponentials, 

= riiliexpVi, where Vi are eigenvectors under the adjoint action of the elements £±, 
namely [i^±, Vj] = uj^'Vi. Note that different orders of the exponentials in the product, in 
general, can give different soliton solutions. These two assumptions concerning a choice of the 
fixed mapping n^, clearly are equivalent to imposing some initial conditions on ( ^.3| ), which 
are relevant for obtaining iV-soliton solutions of our system. The solutions are characterised 

by soliton parameters u;^*^ = cui. — uj_ and f = — r-^ — r-r, and, if one interpretes v^^^ as 

Uj)^' —UJ_' 

velocity of the i-th soliton, cj^-'tui*'* must be negative, otherwise |f *-*''| > 1; moreover, velocities 
v^^^ can be taken to be equal one to each other. To be more precise, 

N 

Tmnifi') = ^'"^ < J]e^»^^P[-'''(^-^'''*)l|/l (4.4) 
1=1 

where we are back to the spatial and time variables, x± = t ± x. The final step in the 
construction of the soliton solutions is to use an appropriate basis vectors \ h >^'^^ and < h'\ 
in the representation space, annihilated by all the subspaces Q+s and Q-t, s > n, t > m, 
respectively. 



5 Dressing transformations 

Now we present an alternative way of constructing solutions given by the solitonic speciali- 
sation of the general solution discussed above. They are obtained by the dressing transfor- 
mations which are non local gauge transformations acting on the gauge potentials and 
leaving their grading structure invariant j^, |10|, [12|. The dressing procedure requires ex- 



istence of two gauge tranformations mapping a given to the potential with the same 
grading spectrum. Namely, there must exist G+ and 0_ such that 

A^ Ap^ = 0±A^0±i - d,e^eg\ (5.1) 

with Afj_ and A'^ having the same grading structure as ( |2.5| ). Here B± are mappings from Ai 
to G. Since satisfy the zero curvature condition, they are of the form 

A^ = -d^TT-\ (5.2) 



'"^By convention, here and in what foUows derivatives act only on the first term of products, unless 
parentheses indicate otherwise. 
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Equivalence of these two transformations implies that 



TP = e+T = e^Tp, or, equivalently, 01^6+ = TpT'^ (5.3) 

with p being a constant mapping from A4 to G, and = —d^TP{T'')~^. 

The dressing transformations are defined, in fact, in terms of a modified Gauss decom- 
position of the element TpT"^, 

TpT-' = {TpT~^)^^ {TpT^'), (^P^'')>o (^-^^ 

with {TpT~^)^, iTpT^^)^Q and (TpT^^)^^ belonging to the subgroups of G, whose Lie 
algebras, defined in (|2.1|) , are Qq, and respectively. Then, we define the element 

as 

TP = ef (TpT'^)^^ T = 0L°^ (TpT-^y^^ Tp, (5.5) 
where and G^'' are introduced in such a way that 

e^^'ef = (TpT-i) . (5.6) 



Comparing with (5^), we write the elements Q± in (|5.1| ) as 

e+ = efe>, e_ = eL°^ef; (5.7) 

with 

e> = (rpT-i)^^, Q< = {TpT-'y^\ (5.8) 

Therefore, the dressing transformations provide a mapping between solutions. Hence, 
the space of solutions to the system can be splitted into orbits of such transformations. The 
knowledge of a given simple solution is then enough to generate a whole orbit of solutions. In 
the case of the models under consideration, interesting solutions, namely solitons, fortunately 
are on the orbit of a vacuum solution. That is the fact we will explore; it is related to the 
solitonic specialisation discussed above. 

Perform now the dressing transformation by taking as an initial configuration a vacuum 
solution of ( 2.13|) - (|2.16|) . As we have said, the model under consideration may have several 
type of vacuum solutions. However, here we will deal with the solutions of type ( |2.18| ) or 

(PO). 

For the vacuum solutions ( |2.2q ), the gauge potentials ( p.5| ) become 

= -£+ , =S- + Qx+C, (5.9) 

with S± given by ( |4.1| ). They can be written as 

A^^^ = -d±ToT,-\ (5.10) 

with 

To = e^+^+e"^-^-. (5.11) 
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The gauge potentials for the vacuum solution ( |2.18| ) are obtained from ( |5.9| ) by taking 
= 0. In fact, they are connected by the gauge transformation 

= foAf f,' - d±fof,-\ with To = exp[x+ (^+ - El)] exp[-x_ {S^ - E_i)]. 



However, in general, the vacuum solutions (|2.18| ) and (p.20|) may not be connected by any 



dressing transformation, and, in such a case, the existence of two elements of form ( |5.7| ), is 
not always possible. Consequently, one can have soliton solutions lying on different orbits 
under the dressing transformations. 

In order to perform the dressing procedure, we take (|5.9| ) as initial gauge potentials. 
Then, we obtain, under the dressing procedure, the solutions on the orbit of vacuum ( 2.2CI| ), 



and for = those on the orbit of the vacuum ( |2.18| ). From the structure of the dressing 



transformations and from the fact that the grading operator (|2.2| ) is never the result of any 
commutation, since it contains D, it follows that the dressing transformations do not excite 
the field t]. Therefore, from (|3D, ( ^JOD , Q and we get 

bEib-'+Y^b F+ r 1 = e± El + E c+ e± e„ e^^ + 9+e± qi\ (5.12) 

m=l n=l 

l-l 



-d_ bb-^ -d_ {u + nx+x_) C + E_i+Y.Fm 

m=l 

l-l 

= e± E^i + E 0± E^m - 9_e± e^^ (5.13) 



m=l 



Note that in the above relations, the fields b, v and F^ stand for the solutions on the orbit 
of the vacuum solution ( |2.20| ). The procedure to construct the solution requires to split the 



above equations into the eigensubspaces of the grading operator ( [2.21 ). It is convenient to 
write 



0> = exp (E^^'M ' 0- = exp (E^^^'M ' where t^^'^ G 

\s>0 / \s>0 / 



(5.14) 



The mappings t'-^*-' , for each choice of p, are determined from (|5.8| ) with T being Tq given in 
( p.ll| ). Then, the components of ( |5.12| ) and ( |5.13| ) in each eigensubspace, give an equation 
connecting the fields with t^'^^^ . Thus the solutions for the fields 6, v and F^ are determined 
from t'-^^-'. Such a procedure is rather cumbersome, but fortunately, one needs to know very 
few t'^^^'^'s to get the solution. For instance, taking relations ( ^.12| ) and (^.13|) for 6+ (©-) 
with grade components and — / (/ and 0), one gets 

ef = hrj}{xJ) , eL°^ = 6e("+^"+"-)^/iR(x+), (5.15) 

with hi^{xS) and h^ixj^ defined in ( |2.32| ). 
From (|^3D, (|5TTD and ( gTTsD it follows that 



ef~' (/iL(x_)6e("+^"+"-)^/iij(x+)) ' e> = e"+^+e-"-^-pe"-^-e-"+^+. (5.16) 
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The space-time dependence of the r.h.s. of the above relation is given exphcitly. One can 
extract the solutions out of (|5.16|) by taking the expectation value of its both sides between 
suitable states of a given representation of ^, in a similar way to that one explained in section 

B 

The solitons solutions are obtained from ( ^.16| ) by choosing the fixed group element p, 
characterising the dressing transformation, as the exponential of an eigenvector of £±, i.e. 

p = e^. (5.17) 

That is the solitonic specialization discussed in section ^. Indeed, if V satisfies the relations 

[£±,V]=uj±V, (5.18) 

then (|5.16|) reads as 

exp (^e^'+'^+-^'-'^- v) = exp (^e^(^-''*) v) , (5.19) 

with 7 = a;+ + a;_, and v = (a;_ — a;+) / (a;+ + since x± = t ± x. 

Therefore, for each eigenvector V, expression (|5.19|) corresponds to a solution that travels 
with a constant velocity v without dispersion. Depending upon the properties of V, as we 
will see below in the examples, such solutions correspond to one-soliton solutions. 

The multi-soliton solutions are obtained by taking p to be the product of several one- 
soliton p's, i.e., 

p = e^ie^2e^=' ...e^^, (5.20) 

with each Vi satisfying [£"-1- , Vi\ = u]^Vi. 

Notice that, under the global gauge transformations ( p^.33D , the gauge potentials 



are transformed as A± hr, A± hjj^. Therefore, from ( ^.2|) one has T — * huT, and conse- 
quently ( 1^ ) implies — > hDQ^fi^ and 9f hjjQ^K^ . Hence, with solution ( p. 161 ) 



corresponding to a fixed element p, a solution, obtained from that by a global gauge trans- 
formation (|2.33|) , is given by ( |5.16|) with the replacement 



p^hnph-o\ (5.21) 

if the condition h£,8± = S± is satisfied. For the solutions on the orbit of the vacuum 
( p.l8| ), that is indeed true, since £± = E±i; see (|4.1| ). For the solitonic case, one then obtains 
for each eigenvector V oi S±, an orbit of equivalent one-soliton (or multi-soliton) solutions 
generated by hijV h]^. 



6 Stress tensor and Hamiltonian reduction 

The two-loop WZNW model was introduced in ; it leads, under the Hamiltonian reduction 
procedure, to the conformally affine abelian Toda systems. The structure of the two-loop 
WZNW model was further studied in ||2^; in the Hamiltonian reduction was extended to 
non abelian affine Toda models. The procedure we discuss here is, in fact, simpler than those 
in since, due to the extra higher grade fields of the reduced model, all the constraints 
are of the first class. 
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We use the Hamiltonian reduction to obtain, from the Sugawara construction for the 
two-loop WZNW model, the energy-momentum tensor for the systems defined in section 
0. For such systems it is not easy to get the canonical energy-momentum tensor in a direct 
way, because, in general, the expression for the Lagrangian is not known yet; and, in fact, 
in terms of the fields introduced in ( p. 61) , it is non local. 

The action for the two-loop WZNW model is the same as that for the ordinary WZNW 
theory, with the fields g being mappings from Ai to G. Therefore, the equations of motion 
for the two-loop WZNW model are 

d+ (d^gg~')=0 , a_(r'5+^)=0. (6.1) 

We are interested in those mappings which can be represented in the form 

g = NBM, (6.2) 

where N, B and M are generated by the subalgebras Qq and respectively, defined 
in (|2.1|) . Introduce the mappings K^in as 

d^gg-^ = N KlN-^ = N (^N~^d^N + d^B B-^ + Bd^MM-^B-^^ N-\ 
g-^d+g = M-^KrM = M-^ [B-^N^^d+NB + B-^d+B + d+M M-^)m, (6.3) 

and so, from (|6.1[), one obtains 

d.Kn = -[KR,d-MM-^l 

d+KL = [Kl, N-'d+N]; (6.4) 

cf. those in p, H. 

We show now that the system ( ^.Tj ) - ( p. 91 ) can be obtained from the Hamiltonian reduc- 
tion of the two-loop WZNW model by imposing the constraints 



(g-'d+g)^ = El, [g-'d+g)^^ = 0. (6.5) 

From ( |6.3| ) one sees that constraints ( |6.5| ) are equivalent to the following ones: 

(9_MM-^)_^ = B'^E^iB , (9„MM-^)^_^ = 0, 
(^N-^O+n)^ = BEiB-\ [N-^d+N)^^ = 0. (6.6) 

Therefore, the mappings d-MM~^ and N~^d+N have components in the subspaces Q-n 
and On, respectively, with 1 < n < / — 1. To provide a relation to system ( |2.7|) - (|2.9|) , denote 



i-i «-i 
Bd-MM-^ B-^ = E.I+ J2 F- , B'^ N'^d+N B = Ei + J2 K = ■ (6-7) 

m=l m=l 
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Substituing the mappings Kl/h, defined in (|6.3|), into ( |6.4| ), and using constraints 
one can easily check that B, and F~ satisfy (p^)-(pl9|). Hence, the system can indeed 
be obtained by the Hamiltonian reduction from the two-loop WZNW model. 

Similarly to the finite dimensional case of the standard Toda system and its higher 
grading generalizations Q, one can prove that the mappings A^, M and B entering the 
decomposition g = NBM, are determined by the holomorphic mappings fi±, see ( p.3|) , via 
the formulas 

N = a+(x_)(7o-)-V+7o- , M = a_{x+)i^+)-'u_^+ , (6.8) 

and the equality 

So, the WZNW field g is represented in a holomorphic factorisable form 

= 9L{x+)gR{x.), (6.9) 

with 

gR = fi^{x^)%{x-)a-\x-), gl^ = fi+{x+)-f^{x+)aZ\x+). (6.10) 

Here 7^ are arbitrary mappings determining the general solution to our system; a± are 
additional arbitrary mappings generated by the subspaces J2m>o G±m, specific for the WZNW 
theory, and they become absent (or, in a sense, hidden) after the Hamiltonian reduction to 
the Toda type theories. 

It follows from equations ( |6.1| ) that the corresponding chiral currents, 

Jr{X;x+) = -k Tr (Xg^'d+g) , Jl(X;x_) = A; Tr {xd^gg~') , (6.11) 

with X E Q, satisfy the equations d^jR = d^Ji = 0. 

The two-loop WZNW model is conformally invariant, and its energy-momentum tensor 
is given by the Sugawara construction [16, 17|. On the classical level, we have its components 
given by 

T++ = ^ Tr d+g)\ T__ = ^ Tr (9_ g g-') ' , (6-12) 
and T+_ = T+_ = 0. Here k is the central term of the two-loop current algebra, 

[J(X, x) , J(X', y)U = J([X , X% x)5{x - y) + kTi {X X') 5'{x - y) ; (6.13) 

X and X' are elements of an affine Kac-Moody algebra Q. The above relation is satisfied by 
the left and right currents (|6.11| ); and currents of different chiralities commute. 

The components T++ and T have vanishing Poisson bracket denoted [*, *]pb, and each 

of them generates a copy of the centerless Virasoro algebra, 

[T(x) , T(y)]pb = 2T{y)5\x - y) - T\y)5{x - y) . (6.14) 

The left (right) currents have vanishing Poisson bracket with (T ), and are trans- 
formed under T (r++) as primary fields of conformal weight 1, 

[T{x) , J(X; y)]pb = J(X; y)5\x - y) - J'(X; y)5{x - y) . (6.15) 
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Under constraints ( |6.5| ), the currents Jr{X^ and Jl{X^'-^; x), such that 

Tr (x(-')Ei) 0, and Tr [x^^^E_i) 0, take fixed non vanishing values. Therefore, the 
Virasoro generators do not weakly commute with such currents. In other words, since the 

currents are not scalars under the conformal transformations generated by T++ and T , 

it means that constraints (|6.5|) break the conformal invariance. However, one can use the 
currents in the direction of the grading operator Qs given in (|2.2| ), to improve the Virasoro 
generators. Define 

L++{x+) = T++{x+) + jJ'RiQ,,x+) , L__ix_) = T__ix_) -jJ'l{Q,,x_). (6.16) 

One can easily verify that these components generate two commuting copies of the Virasoro 
algebra, which are centerless because Tr((5s)^ = 0, for the particular grading operator in 
( p.2|) . In addition, these generators weakly commute with constraints (|6.5|) , and, therefore, 
the reduced model is conformally invariant. Substituting constraints ( |6.5| ) in ( |6.16| ), we 
obtain the energy-momentum tensor for the reduced model, which generate two commuting 
copies of the Virasoro algebra under the Dirac bracket. Indeed, one can easily check, using 
(PD, ( |6A2D , ( leTID and (|6j), that 



ii:Tr((9+F+)„(M(3.M-l„ 



n=l 



+ E (l - y) Tr [{m-'O+m)^^ (m-'F+m)J , (6.17) 



and 



^ n=l ^ 

+ E(i-T)Ti'((a-iViv-'),.(A'f-A'-)J. (6.18) 





]V = exp VC. , Af = exp y;^_. , (6,19) 



Representing 



one sees that the mappings F"*" and F~ defined in ( p.TD , depend on (g and respectively, 
for 1 < s < / — 1. In addition, due to the grading structure of the terms entering ( |6.17D 
and ( |6.18|) , the fields (s and C,-s, for s > I, do not contribute to U^'^ and Uf^. Therefore, 



the components ( |6.17|) and ( |6.18| ) of the energy-momentum tensor, are local functions of the 
fields B, (s and 1 < s < / — 1, of the reduced model. If one wishes to express those 
components in terms of F^, one gets a non local expression. That is, in fact, a difficulty for 
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obtaining the canonical energy-momentum tensor for system ( |^.7| )-( P^ . The corresponding 
Lagrangian, written in terms of F^, is also non local. 

The canonical energy-momentum tensor ( |6.17[ ) - ( |6.18| ) is conserved and traceless as a 
consequence of the fact that it is the reduced two-loop WZNW stress tensor, i.e., 



0, 



red. 



0, 



' red. 



^r_cd. 



0. 



(6.20) 



Notice that such a tensor has a part which is a total derivative, namely 
Tr (Qs d+ (B-^d+B)) in ( |eT7|) , and Tr (Q^ 9_ {d^B B-^)) in They do not correspond 

to the full reduced improvement terms J'n/iiQs), but only to part of it. We can use them to 
construct a trivially conserved tensor, namelyQ 



S^, = --Tr (g, [d^ [B-'d, b) - g^,d, [B'^d^ B 



I 



{d^d,-g^,d^)u, (6.21) 



where, in the last equality, we have made use of the fact that, thanks to the basis chosen to 
parametrise -B, see ( |2.1(j| ) - ( |2.12D , Qs is orthogonal to all generators of B, except the central 
term C. So, such a tensor is symmetric and conserved]] 



d" S, 



0. 



Then, introduce an energy-momentum tensor B^jy as 



e 



T rod. c 



Due to (|6.2CI| ) and ( |6.22|) , it is also symmetric and conserved. 



9^6 



fii/ 



(6.22) 



(6.23) 



(6.24) 



but it is not traceless. 



7 Masses of fundamental particles and solitons 

As we have seen above, the system under consideration is conformally invariant. There- 
fore, since we do not have a continuum mass spectrum, its fundamental particles have to be 
massless. However, such a symmetry can be spontaneoulsy broken by choosing a particular 
constant solution for the field r], say t] = t]q. The resulting theory is then massive. Repre- 
senting the mapping B as B = exp T, and considering only the linear field approximation, 
i.e., the free part of the equations of motion (p.7D-(py9D, one gets 

d+d^T = -Vr,[E.i, [Ei,T]], (7.1) 
= -v,[E,i , [E, , F^]] , (7.2) 
d+d^F- = -v,[E.i , [E, , F-]] , (7.3) 

^The metric is g_|-+ = g = 0, = g ^ = 1/2, and since B commutes with Qs, one has 

9^Tr (QsB-^d^B) ^ Tr {Q^d^BB-^). 

^Due to the fact that B commutes with Qs, it foUows that S*^^ is symmetric and conserved independently 
of the basis we use. 
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where = e^^°. 

Therefore, the masses of fundamental particles in such a theory are given by the eigen- 
values of the operator [E_i , [Ei , *]] in the subspaces Qn, n = 0, ±1, ±2, . . . ± (Z — 1), i.e., 

[E_i, [Ei,X]] = XX. (7.4) 

Since = ^{df — d"^), we obtain the masses from the Klein-Gordon type equations (|7.1|) 



- (|7J) as 

ml = 4:Xvrj. (7.5) 

That result constitute a generalization of the arguments used in the abelian and non abelian 
affine Toda models |^. Of course, we are interested in those cases where the eigenvalues of 
the operator [E^i , [Ei , *]] are real and positive on the subspaces under consideration. That 
will be, in fact, one of the conditions we use to select the data {Q,Qs,l,E±i} for defining 



physical models through (2.5) 



Notice that the field plays the role of a Higgs field, since it not only spontaneously 
breaks the conformal symmetry, but also because its vacuum expectation value sets the mass 
scale of the theory. We have here the same mechanism as in non abelian affine Toda theories 



The masses of the fields (m and C,-m, introduced in ( |6.19| ) for 1 < m < / — 1, are the 
same as those of the fields F^. Indeed, substituting ( |6l7|) in ( |2.14|) -( p.l5|) , and taking the 
linear field approximation with t] = t]q, we get 

d+d^Cm = Vn[Ei, [E^i, Cm]], 

d+d-^-m = V^[Ei , [E.i , ^.m]], (7.6) 

where a trivial integration is performed to eliminate one derivative. 



Let us explain now, following the reasonings of |^ and g , that the masses of solitons are 
also generated by the spontaneous breakdown of the conformal symmetry. The energy of clas- 
sical solutions are given by the space integral of the (0, 0) component of energy-momentum 
tensor U^^- defined in (|6.17|) - (|6.18|) . In the Lorentz frame where the classical soliton solution 
is static, the energy should be interpreted as the mass of the soliton. However, since the 
theory is conformally invariant, it has no mass scale, and the soliton mass should vanish. 
When the conformal symmetry is spontaneously broken by choosing a particular constant 
solution for the field 1], we obtain a massive theory. Construct the energy-momentum tensor 
of such a theory as follows. Clearly, the tensor 9^,^, introduced in ( |6.23| ) and evaluated at 
any classical solution, satisfies ( |6.24]) . Therefore, the tensor defined by 

pvbrokcn _ i (7 7\ 

^ — ^fii/ |r;=constant ) V / 

is symmetric and conserved, 

^/.Qbroken _ q ^ (7 

since rj = constant is a solution of the equations of motion. Then, let the energy in the 
massive theory be proportional to the space integral of Gjj™'^^'^. Using ( |b.21| ) and ( |6.23| ), we 
obtain the soliton mass in the form 

M / \ k N 

' j JxQ'^t"''-E..c.)= j^d^{u + nx+x^) (7.9) 
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because the integral of Uq^' vanishes by the above arguments. Here v is the sohton velocity 
in the units of the speed of the light. Notice that we have subtracted the energy i^vac. of the 
vacuum solution which is, in fact, divergent. Of course, the vacuum solution is not unique, 
and it is not clear which one provides the absolute minimum of the energy. We will use the 
following prescription for the soliton mass formula. For the soliton solutions lying, under 
the dressing transformations, on the orbit of the vacuum solution (|2.20| ), we take Q in ( |7.9| ) 
to be that one given in ( p.21| ). However, for those soliton solutions lying on the orbit of the 
vacuum ( |2.18| ), we take Q in ( |7.9| ) to be equal to the parameter /3 introduced in ( |2.17| ). Such 
a prescription guarantees the finiteness of the soliton masses. 

The soliton masses are determined solely by the behaviour at x = ±00 of the space 
derivative of the field u. That is quite a remarkable fact. In addition, as we now explain, it 
is very easy to obtain such a behaviour in the general case from the solitonic solutions 
or ( ^191) . 



Consider the integrable highest weight representation of Q with highest weight state | Ag), 
satisfying the relations |]T5|, ^ 




\\s) = sa |As), C |As) = ^ E^f^O l^s), f^ \ K) = , if = 0; (7.10) 



cf. (|3.14|), where fa = E^^ , a = 1,2, . . .r, fo = E^^ , and Sj are the co-prime integers labeling 



a given integral representation of Q (^3), if are the integers in the expansion = J2a=i 



and /q = 1. Such a highest weight state can be realised by 

I As)=(g) I A,)®^% (7.11) 

1=0 

where | Aj) are the highest weight states of the fundamental representations of Q, and Aj are 
the corresponding fundamental weights of Q. 

Consider now an integral gradation of Q, with s[ = ^ Si, ao = —tjj, and Si labeling 

the gradation that defines the model ( |2.13| )- (|2.16|) . Therefore, it follows that the Cartan 



generators of the special basis introduced in ( p.llD -( p.l2D provide 

H'a I A.0 = 0. (7.12) 

From ( [7.10|) one has that | Ag') is annihilated by all negative root step operators with s'-zero 
grade, and, consequently, with s-zero grade too. Since it is a highest weight state, | Ag') is 
annihilated by all generators of the subgroup Gq. Therefore, taking the expectation value of 
both sides of ( |5.16| ) in such state, one gets (with the gauge choice hi^X-) = hR^x^) = 1) 



iu+nx+x.)N,^ _ ^^^^ I ^x+£+ ^-x.e- p^x.e- ^-x+e+ | ^^^y ^^^^g^ 



e 



Now, choosing p to be the exponential of an eigenvector of £±, 

[£± , V] = UiV, (7.14) 
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we obtain a soliton solution 

(As' I e^^ ^ I As') (7.15) 



with r = a;+x+ — a;_a;_ = 7 (x — ft). 

Suppose V is an operator in such a representation for wich there is a positive integer Ny, 
such that 

(As' I V I As') =0 for n > N{,. (7.16) 



Then the sohton mass is easily obtained from (|7.9|) , where for 7 > (7 < 0) only the upper 



(lower) limit x = 00 (x = — 00) contributes in the integral^ , 

Notice that we must have uj+uj^ > in order to have the soliton velocity v = (tU- — 
u;+)/(ct;„ + c<j+), not exceeding the light velocity (c = 1). 

The soliton mass formula (|7.17|) has some remarkable properties. One of them concerns 
the relation particle-soliton in the theory, indicating some sort of duality similar to the 
electromagnetic duality of some four dimensional gauge theories possessing the Bogomolny 
(monopole) limit [^]. As we have seen, the soliton solutions are created by the eigenvectors 



V of £^-|-. From ( [7.14|) one has , [S^ 1 ^]] = ^^+^^-V. Expanding V over the eigenvectors 
of the grading operator QsasV = Y,n V^"'\ one observes that , , V*-"-*]] = lo^lo^ 
Therefore, if some V^""^ E Qn, n = 0, ±1, ±2, . . .± (/ — 1), does not vanish, it implies that V^"^ 
must be one of the eigenvectors X in ( [7.4|) . Then we associate a soliton with a fundamental 
particle. In addition, we have A = u+U-, and, consequently, from (fTsP and ( [7.17] ), the 



masses of the corresponding soliton and fundamental particle are determined by the same 
eigenvalue. In fact, we have from ( [7.5| ) and ( [7.17| ), with = 1, that 

Mso, = ^^mr'-. (7.18) 

Of course, in the expansion of V, we may have more than one non vanishing V^'^\ with 
n = 0, ±1, ±2, . . . ± (/ — 1). Then we would associate a one-soliton solution to more than 
one fundamental particle. The counting of one-soliton solutions has to be better analysed 
in each particular case. We discuss this issue in the section devoted to the examples. 



8 Physical properties of the higher grading fields 

It is clear from ( [^.23| )-( p.27| ), that the massive fields associated with non vanishing grade 
(namely F^), are chiral fields with non vanishing spins, in contrast with the Toda type 
fields. In fact, we show that the free equations for such fields take the form of the massive 
Dirac equation, as could be expected from general covariance arguments. 



^We point out that the sohton mass formula (7.17) could be equally obtained by defining the mass through 
the momentum formula, instead through the energy like in (7.£), as Jf^^-^ = / c?x8q™'^°". In this case, we 



do not have to subtract the vacuum momentum, since it vanishes. 
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Consider the subspace Qm for < m < /. Let Qj^^ be the subspace of Q^, generated by 
the eigenvectors of [E_i , [Ei , •]] with non zero eigenvalues, i.e., 

where A*^™'' is defined as 

[E^i , , T^"^)]] = , [E_i , T^'")]] = A^'") T^™). (8.2) 
Decompose the subspace Qrn, as a vector space, into the sum 

g^^g^^^ + g^^\ (8.3) 

where Qj^^ is the complement of in Define now 

Since A^"') ^ 0, it follows that T^-i+"^) is necessarily non vanishing. Therefore, we have 

[El , [E_i , T(-'+-)]] = [E_i , [i?, , T(-'+™)]] 

= [E_i , , [E_i , r("^)]]] = A^'") [E_i , r(™)] = A^"*) r(-^+'")(8.5) 

So, whenever we have an eigenvector in with eigenvalue A^"*^ 0, we also have a corre- 
sponding eigenvector in Q-i+m with the same eigenvalue. Notice that if A'-™-' is degenerate, 
then the corresponding eigenvectors are mapped bijectively. Suppose T^"* and T2™'' have 
the same eigenvalue, and that they are mapped onto the same element in ^_/+m, i-e., 

T^^] = Tf)] , and so T^^ - T^^] = 0. (8.6) 

However, this relation is in contradiction with the following ones: 

[El , [E_i , T^^ - ri"^)]] = A("^) {rt^ - Tt^) , and t'T^ - T^"^ ^ ; A^"^) ^ 0. (8.7) 

Analogously, we write 

G-l^m = G^?+m + G^-ilm, (8-8) 

with 

g^Zm = {T^-'^""^ e g.l+rn \ A^^^ ^ 0}, (8.9) 

and g^_^l„i being the complement of ^i^|„j in g^i^„i- Therefore, from the considerations 
given above, we conclude that the subspaces ^i^^^ ^nf "* isomorphic. The mapping 
is given by the action of E_i on Our arguments are applied equally well in the reversed 

^ ( F) 

direction, the mapping can also be given by the action of Ei on g_i\_^; and hence 

C^] = gZi^, [^-/ , C^] n ^ifi^ = {0}, (8.10) 
[Ei^g^V-^j = C^ [^^ , ^5iJ n = {0}. (8.11) 
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One has 

[E., , Gl^^] n g^Zn. = m , [El , ^ifiJ n = {0}, (8.12) 

because, otherwise, it would be in a contradiction with the relations 

[El , [E.i , g(^^]] = , [El, [E^i , ^iflj] = 0. (8.13) 

Notice that Qm (^-/+m) does not lie necessarily in the kernel of E^i (Ei). However, if a given 
element J^""^ of is in the image of Ei; then it follows from ( p.l3| ) that it must be in the 
kernel of E_i, i.e., 

J^'") = [El , =^ [E^i , J^"^)] = , J(") G (8.14) 

Analogously, 

= , J^"^)] ^ [E; , = , J(") e gi^\ (8.15) 

In the same way as in section |^, we consider the linear approximation where t] = t]q, b = 1, 
u = (recall equation (|2.1CI| )). This gives the free part of equations (|2.8| ) and (|2.9| ), 

= e('-)''°[^/ , Fr_J, d^Er_^ = -e^^-[E_i , (8.16) 

Denote the generators of and ^i^l^; corresponding to the eigenvalue A*^™^ 7^ 0, as T^™^) . 

and T^^^\ respectively, where the index i stands for a possible degeneracy of A^"^). The 
basis is chosen in such a way that 

[E.i , T^S),,] = r^ri^r^ [El , T^rJ^T^)] = TaM tS,. (8.17) 

We also use the notations 

^m—^m ^^ra ' ^ l-m — ^ l-m ^ ^ l-m i 1,0. iOJ 

with 

F+,(F/K) n{F/K) p-,{F/K) MF/K)_ ,^ 

and 

Since the action of E±i does not mix the subspaces of indices F and K, we can split equations 

( prr^ ) to get 

Introduce 

^ = I J , (8.22) 
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and 



7o = 



satisfying anticommutation relations 



{7m ' ^u} 



7i = 




(8.23) 



(8.24) 



with (7oo = —(711 = 1, goi = 0. Therefore, equations ( |8.21| ) can be written as the Dirac type 
equations 

where, following (|7.5|) , we put 



1 + 75 



+ 



75 ^(l-m)rjo 



(8.25) 



(8.26) 



Hence, the massive degrees of freedom of the model, corresponding to generators of grade 
±m with < m < I, are, in fact, Dirac fields. However, in general, the mass term involves 
a 75 term, so that parity in violated. A noticeable exception is when / = 2, m = 1. 

Notice that "ip^/j^'^ are the projections of ip^'^'^^'^ under ^i^^, with 75 = 70 71. Under the 

gauge transformations ( p.30|) -( p.31j ), ip^ and ip^ are transformed independently under 
the action of (2;+) and hi (a;_), respectively. 

At this point we have seen that the mapping between Qm and Q-i+m displayed at the 
beginning of the section precisely ensures the existence of both chirality components of the 
Dirac fields. The next question is whether we may write a free action. Suppressing indices 
for a while, a free field actions would take the form 



Writing 



we obtain 




^ = (^«, ^l)7o, 



C/i 



R(- 



(8.27) 



(8.28) 



^.29) 



where niR = a + b, itll = a — b. Next, equation ( p.25| ) shows that the masses vanish in the 
limit r/o — > cxD. Then, we may discuss the conformal properties of the ip fields in the usual 
language, where a field A{z^,z_) is primary with weights A+,A_ if A{dz^)^+ {dz_)'^- is 
invariant by conformal transformations. The kinetic term of equation ( ^.29|) is conformal 
invariant if we have the following weight assignements. 



A+ 

A_ Jl 1-Ji 



^.30) 
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Consider next the mass term. With the Lorentz transformation x± 6'^^x±, 



and ipi'ipR will be invariant if Jl + Jr = 1- Returning to our Dirac fields, one finds the 
identification 

for some choice of j, A;,£. Thus the fields of type R (L) are to be split between ip and ip 
type fields. The spectrum of weights is just right to form the needed quadruplets provided 
the number of R (L) fields is even. This is not always true as we shall see. Another difficulty 
is that in general if) is not the complex conjugate of ip. Thus the above free action is not real. 
This was already the case for the usual affine Toda theories beyond sine-Gordon, and we 
may expect that this will not be a problem. Another concern is the statistics. At first sight 
our ip fields are c-number fields so that they seem to describe bosons. However, it is well 
known that in two dimensions the statistics of fields depends upon the coupling constant. 
Perhaps the latter should be fixed so that the tp fields become anticommuting operators. We 
will have more to say on this below. 



9 A special class of models 

We now describe a class of models possessing a U{1) Noether current, which, under some 
circunstances, is proportional to a topological current. That occurs for those models where 
the grade / of the operator Ei, introduced in (|2.6| ), is equal to the integer A^s defined in (|2.2|) . 
So, throughout this section we have I = Ns. 

The calculations become simpler if we realize the generators of the affine Kac-Moody 
algebra Q, in terms of those of the finite simple Lie algebra Q as 

H: = z^Ha, E2 = z^E^, D = z^ (9.1) 

with z being a complex variable. Then, the Lie bracket and trace form on Q are given by 

[A{z) , B{z)] = [A{z) , B{z)]g + C f^jT (^B{z)^^A{z)^ (9.2) 

and 

Tr {A{z) B{z)) = / -^tr (Aiz) Biz)) (9.3) 
J zmz 

where [■ , -Jg and tr (■ ■) are the Lie bracket and trace form, respectively, on Q. 

As a first step, we construct two chiral currents associated to the elements z'^Ejq^ and 
zE-N^ of Qo- Multiply eq. (|2.8| ) by z~^E^^_^, sum over m and take the trace. Then, take 
the same equation with m replaced by Ng — m, multiply by z~^F^, sum over m and take 
the trace. Add both to get 

Y: 5_ Tr _ J = -2 5: Tr {z-'E^^[F+ , B-'F-B]) + X+ (9.4) 



m=l m=l 
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where 

Ns-2 Ns-m-l 

= E E tt{z-'fi_jf;1^^,b-'f-b] 

m=l n=l 
Ns-1 m-1 

+ ^ T.T^^{^''K[F'l+n-m,B-'F-B]) (9.5) 

m=2 n=l 

Notice, the first sum in m ends at A^s — 2 and the second starts at 2, because the quadratic 
term in eq. ( |2.8| ) does not exist for m = Ng — 1. 

Similarly, multiply (|2.9| ) by zF^^_^, sum over m, and take the trace. Then, take the 
same equation with m replaced by A^s — fn, multiply by , sum over m and take the trace. 
Add them to get 

E 9+Tr(^F^F^^_^) = 2 ^ Tt{zE_n^[F;^,BF+B-'])-X- (9.6) 

■m=l m=l 

where 

E E 

m=l n=l 

Afs-1 m-1 



^" = E E Tr{zF^^_jF-^^,BF^B-'] 

m=l n=l 

Ns~l m-1 

+ E E Tr {zF- , (9.7) 



m=2 n=l 



Using the fact that 

Ns~2Ns-m-l Afs-lm'-l 

E E ~ E E ' where m' = m + n (9.8) 

m=l n=l m'=2 n=l 

one gets[| 

X+ = X- = (9.9) 

Projecting (|2.7|) onto z~^En^ and then onto zE^n^, and comparing with (|9.4|) and (|9.6|) , one 
getf] 

9_J^ = 0, d+J = (9.10) 

where 

JVs-l 

2 



Ns-m 



1 iVs — i 

:r(xH.) = Tr (z-'E^B-'d^B) - - E Tr (^-'F+F,^ 

m=l 
1 Ns-1 

J{x^) = TT{zE_^^d^BB-')~-Y.TT{zF^F^^_Ji (9.11) 



2 m=l 



^In fact, each double sum in X"^ vanishes separately. Use that J2ni=2 ^"i=i ~ J2n=i'^ ^m=n+i^ 
notice that the terms for (fixed n) m ~ n + k and m = Ng — k cancel. 

^°In the projection onto z^'^Em^, it is easier to use the equivalent form of (2/7), 9_ {^B^^d^B^ = 

-[E^^ , S-i E^N^ B] - Y.n=i'[Fi , B-^ F- B] 



30 



We now consider those models where z ^Ej^^ and zE^^^ are parallel, and lie in the center 
of^o,i.e.H 

zE-N^ = fj.z~^Eiy^ = Eq e center of Qo (9.12) 



where fx is some constant independent of z. Using such condition one observes, from (|9] 
and (|9.6| ), that the current^ 



^+ = E Tr {z-'f:^FI_Ji , J- = - E Tr {zF-F^^_Ji (9.13) 

m=l m=l 

is conserved 

= (9.14) 
In addition, the condition ( |9.12| ) implies that the current 

J+ = - Tr [Eo B-^d+B) , J_ = Tr (^o d^BB-^) (9.15) 

is a topological current, i.e. it is conserved independently of the equations of motion 

d^r = (9.16) 

We now come to a very interesting property of these models. Under the conformal transfor- 
mations ( p.22|) -( p.27|) , the chiral currents (|9.11|) transform as 

J{x+) (/(x+))-'(^^(x+)--^Tr(EoQs)(ln/'(x+))'j 

J{xJ) ^ {g'{x^)y'{j{xJ)~^Ji{E,Q,){\ng\x^))'^ (9.17) 

Therefore, if 

Ti{E^Q,)^Q (9.18) 

one concludes that, given a solution of the model, one can always map it, under a conformal 
transformation, into a solution where^ 

J{x+) = J{xS) = (9.19) 

Such a procedure amounts to a gauge fixing of the conformal symmetry. We are choosing 
one solution in each orbit of the conformal group. Another way of saying it, is that ( |9.19|) 
are constraints and we are performing a Hamiltonian reduction. The degree of freedom 
eliminated corresponds to the field r]. So, in the submodel defined by ( 9.19| ) one observes 
from (|9.11|) , that the Noether and topological currents (|9.13|) and ( |9.15|) are equal 

J, = 3, (9.20) 



"'^"'^ Obviously, Eq can not have any component in the direction of the central term C or the derivation D, 
since it is the projection of elements belonging to Q±n^- 

^^We use symbols with a tilde since the currents we are defining right now are not yet properly normalised 
(see below). 



13 



The exceptions occur when J{xj^) and/or J{xJ) present singularities. 
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as a result of the field equations. As we will see on explicit examples, the true Noether and 
topological currents differ from the one just defined by overall constants: 

J fj. Cpvjoethor^/^) Ctopol.J/i (9.21) 

Thus the topologican and Noether currents are proportional 

J, = ^J^^j^. (9.22) 

Ctopol. 

This is a very important property of such models, which can lead, in some cases as we will 
discuss below, to the confinement of the matter fields. In general, that is if JT" 7^ 0, or 7^ 0, 
one gets from eqs.( |9.1lD , (|9.13|) and ( |9.15|) , 

I2J = J+-U J = -J-+ 3- (9.23) 
Therefore the Noether and topological charges 

/CxD poo 
dxJt gtopoi. = / dxjt (9.24) 
-00 J —00 

satisfy 

Q'Noether Q'topol. 1/ \ /nr>rN 

— = 7^{(lR- Ql) (9.25) 

^Noether Qopol. ^ 

where we have denoted 

/oo roo _ 

dxj{x+); Ql = dxj{x^) (9.26) 
-oo J —oo 

We now comment on the relationship between the masses of particles and solitons of the 
theory and symmetries associated to Eq. From ( p.l2| ) one gets that Eq commutes with -E±Ar^ 
and therefore it generates a diagonal U{1) global gauge symmetry of the type described in 
( p.33| ), namely 

with 9 = constant. The charges of the fields associated to such U (1) symmetry are, of course, 
given by the eigenvalues of Eq in the subspace ^fidds = 0^=~Ars+i ^n, i-e. 

[Eo,T] = qT, Te ^fields (9-28) 

The masses of the fields are determined by the eigenvalue of [Ej^^ , [-E-at^ , ■]] (see ( |7.5|) ). So, 
using ( p.l2| ) 

[E^^ , [E_N^ , T]] = ^[Eo, [Eo ^T]] = jT (9.29) 

and so 

m| = - (9.30) 
A* 
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The solitons of such theory are created by the operators 



VTiC)= E r"T^ (9.31) 

n=— oo 

since they are eigenvectors of E±n^, (see ( ^.16|) -( ^yT9l) , T" = z"-T) 

[En^,Vt{C)] = -qVTiO (9-32) 
/i 

[E^N., VriO] = C'qVTiO (9.33) 

Therefore, if the expansion of the exponential of Vr(C) truncates, we get, from ( [7.17| ), that 
the one-sohton masses are 

Msoi. ~ ^ (9.34) 

Therefore, the masses of the fundamental particles and solitons of such theory are pro- 
portional to the U{1) charge. We have here a situation similar to four dimensional gauge 
theories with Higgs in the adjoint and in the BPS limit, where the masses of particles and 



monopoles (dyons) are given by mass ~ yQeiect. + ?mag.- That point has to be further inves- 
tigated, because we believe it indicates these systems possess some sort of duality involving 



particles and solitons ||30[| . 



9.1 The example of the principal gradation 

Some special examples of models which possesses the structures described above can be 
constructed as follows. Let Q be any untwisted affine Kac-Moody algebra furnished with 
the principal gradation, where s = (l,l,...,l), and corresponding grading operator, Qppai, 
given by ( |2.2| ) with Ng = h = Coxeter number. Therefore 

Go = {i/°,a = l,2,...r;C;gppai} 

Q-m = (9.35) 

where < m < h, and a*-™-' are positive roots of height m, i.e. they contain m simple roots 
in their expansion. Following ( p.lO| ) we parametrize B as 



where is defined in (p.ll|) , and z/ = z/ — |5 ■ yj, with 5 = Y.a=i ^5 Xa being the 



fundamental weights of G- We then choose E±h to be parallel 

E±h = ni- H^' (9.37) 

where we shall choose m to be a vector inside the Fundamental Weyl chamber (FWC), and 
so m ■ a > for a > 0. Consequently 

EQ = m-H^ (9.38) 
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satisfies ( |9.12| ) (with n = 1), since [Eq , Qq] = 0. Since the masses are determined by the 
eigenvalues of [Eh , [E_h , ■]] (see (|7.5| )) we conclude that the y^'s fields are massless, and the 
fields in the direction of the step operators have masses = 4(m ■ a)^. In addition, since 
m lies inside the FWC, one has rria 7^ for any a. Therefore, according to the discussion 
in section (|^), all the fields in the direction of the step operators are Dirac fields. Then, 
following (|8.17| ) and ( p.2q ), we write 

with < m < h, and where we have denoted 

^^.c™) = m^o.(m) = 2m ■ a^"^) (9.40) 
Consequently, associated to each positive root a^'^\ we have two Dirac fields 

r- . ( ) ; . ( ) (^.d 

Notice such notation is in accordance with (|8.31| ), since V'r'"' ^ind ^ire in the direction 

of -E'°{„i) € Gm and -E'i^(m) G respectively. Similarly, and are in the direction 

of E~l^^ e G-h+m and G respectively. 

These systems possess a ([/(l)^)' ® {U{1)r)^ gauge symmetry of the type ( ^.2!j| )-( P?^ ), 
with 

/il(x_) = e''^("->^" ; hR{x+) = e*^«("+)-^° (9.42) 

since these h^iR satisfy ( p.32|) . They also possess a global gauge symmetry of the type (|2.33| ), 
namely 

^^ip- v^v- v^r^; ^ e^^-^V^e-^-^° (9.43) 

with 6 = constant. Notice that the charges of the fields with respect to U (1) global symmetry 
(^) are (see (|Og| )) 



= % = (!'> = ^ V*"' " = m ■ (9.44) 

The equations of motion for these systems, obtained from ( 2.13 )-( 2.16|) , are 



aV = -4 X: E ^^5(«(-)-^+r.,) /l + Ts _ (ijll^'j (9.45) 

52^ = -2 E E ^%.(™) e'^' e^^(°''"'-^+™'') (1 - 75) V^"""' - 4m2 e'^^ (9.46) 

m=l Q,(m) a*^™^ 
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i^^^d.r'^' = m^,(™)e^^("'"'-^+H (^il±^ + e'^''il_^^^-'"^ (9.47) 
d'^r] = (9.49) 



where the Eramma matrices are defined in (g]^, t/^"^"' = 70, andQ 



T 



u"-'^' = ( ) ; u"-"^' = ( ) (9-50) 



withQ 



(m)2 m— 1 
Y^'^^aC™) n=l 

= E e-Tr(EV)[^„l+n,e--^°F-e-^°]) 

W^^^at™) n=l 
(m)2 h—m—1 

= E e-Tr(E%[F-^„,e-^°i^„+e~-^"]) 

n=l 
„(r?i)2 m-1 

= ^^=Ee"''Tr(E-(L,[F+__„,e--^V-e-^"]) (9.51) 

Y^^^at™) Jl=l 

Since these systems satisfy (|9.12|) , they possess the conserved Noether and topological cur- 
rents ( p.l3| ) and ( |9.15| ), respectively. In order to correctly define these currents, we have to be 
more precise about normalizations. First, concerning the topological current, the potential 
terms in Eqs. |9.45| - b).49 involve only through exp(±a^™').(^). We see that the potential is 
invariant under 

ip^ip + 27ri/i'' (9.52) 

where ii" is a co-weight of i.e. u" = Y1 ^17 ^a'^^ with being the fundamental weights 
of Q, satisfying ^^^2^ = dab, and being the simple roots of Q. Therefore, /i*" ■ a is a 
integer for any root a of Q. We shall choose m = fioJ^a^^a^a, with G 2, and where 
Ho is an overall scale parameter with the dimension of a massj^. Under Eq. |9.52 , we have 
m.ip/ fiQ — > m.ip/ fiQ + 27riZ. The appropriate definition of the topological current is 

f = TT^^^'du (m ■ y;) . (9.53) 



^^Notice that we are not assuming ip" = -0" " 

""^^We use the following normalization for the trace form, Tr {x ■ H^y ■ i?^") = x-y, and Tr (^E^EZ2) = 
^^Notice that ^0 and all Wa's have to have the same sign in order for m to lie in the Fundamental Weyl 
Chamber (see (9.37)) 
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Indeed, the vacua are infinitely degenerate and for a soliton solution the asymtotic values of 
ip, a.t X = ±cx), have to differ by values appearing on the right hand side of Eq. |9.52| . Thus 



it follows from the argument just given that Qtopoi. = / dxj E Z. Concerning the Noether 
current, we have to take account of the fact that, with our field normalization, the free part 
of the Lagrangian has an overall factor k, where k is the coupling constant. The correct 
definition of the Noether current is 



j, = ^ E E ^^.c-.r'-'T.r""' - ^ E ri.r 0.54) 

A^O m=l aim) a^"") fJ'O a>0 " 



where a stands for any positive root of Q. Using the special form of m just introduced, we 
obtain 

2 

= E E ma—a.aa kip^-ff, (9.55) 

«>0 a ^ 

For each term A^" = / kip"''yo is such that 

Moreover, as is well known, ^a.Oa G Z. Thus, the above definition of the Noether current 
is such that the Noether charge ^Noether = / dxJo has integer eigenvalues, as it should. 
Next let us compare the two currents so defined. Clearly Tr (Qppai-Eo) = 25 ■ m 7^ 
=1 tJ")' since m lies inside the FWC. Therefore ( p.l8| ) is satisfied, and one can always 
find one solution in each orbit of the conformal group such that ( p.l9| ) is true. Consequently, 
for such solutions we have 

i E = ^^"d, (m ■ if) (9.56) 

Consequently the Noether and topological currents are proportional: 

3^ = Qtopol. = -^QNoether (9.57) 

This equation has an important consequence. It will certainly hold at the quantum level, 
after a suitable redefinition. Since the eigenvalues of the Noether charge will be integers, 
and, since the topological charge is also an integer, we obtain that kir should be rational. 
One may expect that this will lead to the correct statistics for the ip fields, following the 
argument given at the end of the previous section. A quantification of k is of course expected, 
since our actions are related with the one of WZNW. On the other hand, one may regard 



eqs.9.57 as classical versions of a formulae of the Atiyah Patodi Singer type (see e.g. ref. |31 



for a review). In pratice, Eqs. |9.57| mean that the Noether density is non zero only where 



dip 7^ 0, that is inside the solitons. Thus the ip fields are confined inside the solitons. We 
shall come back to this on the simplest example below. 

The soliton solutions are obtained as follows. The operators diagonalizing the adjoint 
action of E±h, given in (|9.37|) , are 

K (C) = E C'^K (9.58) 

neZ 
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for any root a (positive or not) of Q. Indeed 



\E±K , V; (C)] = C^^ (m • a) (C) (9.59) 

Notice that (C) and V_q, (— C) have the same eigenvalues, and it turns out that the one- 
sohton solutions are obtained from ( |5.16| ), by choosing the constant group element p to be 



an exponentiation of (see discussion below (|5.16|) ) 

(al, C) = < (C) + (-0 (9.60) 

We then have a one-soliton solution (with two parameters a^) for each pair of Dirac fields 
and '?/'". The masses of the solitons and Dirac particle are proportional to m ■ a. 
The construction of the soliton solutions, as well as the physical properties of such models. 



are discussed in detail in section [T^ for the case of s/(2). 



10 Example of the principal gradation for s/(2) 



(1) 



In this section we discuss two examples, for / = 2 and / = 3, associated with the principal 
gradation of the untwisted affine Kac-Moody algebra s/(2)*^^\ Let us denote by if", £"5, D 
and C the Chevalley basis generators of the 5/(2)^^-*. The commutation relations are 

[H'^.H^] = 2mC5m+n,o, (10.1) 

[H'^,El] = ±2E;^+", (10.2) 

= if'-+" + mC5^+„,o, (10.3) 

[D,T'"] = mT™, T"' = H"',E!^; (10.4) 

all other commutation relations are trivial. The grading operator for the principal gradation 
(s=(l,l))is 

Q = -H° + 2D. (10.5) 
2 

Then the eigensubspaces are 

Go = {H',C,Q}; 

g^n = {i^n, nG{Z-0}. (10.6) 
The mapping B is parametrised as 

B = e^"' e'^^e''^ = e^^" e^^'e^^, (10.7) 



where = — | C is the Cartan generator in the special basis introduced in ( ^.111 ), and 



so z/ = z/ — \ip. 
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10.1 Case / = 2 

Consider the case 1 = 2, and choose 

E2=mH\ E_2 = m H-\ (10.8) 
where m is a constant. We then have 

[E2, [E.2,El]]=4m'El. (10.9) 

Therefore, each of the subspaces Q±i has two generators with the same eigenvalue 4m^. 
Following ( p.lTp and ( p. 201 ) we write 

F+ = 2v^ {iJn El + i,RE\) , F{ = (jPl E^^ - ^° ) , (10.10) 
and introduce the Dirac fields 

From (|7.5| ) we obtain the masses of the particles, 

= rriy = rrir, = 0; = 4m; (10.12) 

The equations of motion derived from ( p.l3| )-( ^.16| ), are 

= -4m^V'75e''+^'^^=V', (10.13) 

8^9 = -2m^^(l-75)e''+2^^5V- i^Je^", (10.14) 



1 

-1 
2 
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8^7] = 0, (10.15) 
i'yf'd^^ = m^e''+2'^^^^, (10.16) 
i^d^^ = m^e'>-'^'P^'ij, (10.17) 

where the gamma matrices are defined in ( ^.23|) , and 75 = 7o7i, and ip = ip'^ •jq. Recall that 
d"^ = d"^ — d^, x± = t ± X. The corresponding Lagrangian has the form 

+ z^/^7^9/,^/^ (10.18) 

It is real (for 77 = real constant) if ^jj is the complex conjugate of ijj, and if (/? is pure imaginary. 
This will be true for the soliton solution as we shall see below. 
Notice that such model is invariant under the transformations 

x+^x-] i)R^e^L] i)R^-ei)L] V^^V^; 7] ^7]] v^v (10.19) 

where e = ±1. It should be interprated as the product CP of charge conjugation times 
parity. Parity alone is clearly violated. 
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The generator G Qo commutes with E±2, and, therefore, the gauge symmetry ( p^.29| )- 
of the model is U{1)l ® U{1)r, 

hL{x_) = e^-^"-)^" , hR{x+) = e«+("+)^°. (10.20) 

Since the genereators of U{1)l and ?7(1)_r are the same, we have the global gauge transfoma- 
tions ( |2.33|) generated by = hi = h]^ = e*^^"/^ [6 = const.). The fields are transformed 

as 

V'^e'V i)^e~'^i) u^u, 7] ^7]- (10.21) 

and the corresponding Noether current is 

7/^ = ^7/^^, 5^JA' = 0. (10.22) 

The fields ip and ijj have charges 1 and —1, respectively; and z> and rj have charge zero. 

Let us next see how the general arguments given above concerning Noether and topolog- 
ical charges apply here. The topological current and charges are 



3' 



^t^-d,^, Qtopoi= I dxf, (10.23) 



Indeed, the Lagrangian ( [10.18| ) has infinitely many degenerate vacua due to the invariance 
under (p ^ (p + in. Making use of the field equations, one easily verifies that 



(10.24) 



Combining this equation with the conservation of the vector current tp'j'^ijj, one deduces that 
there exist two charges defined by 

which satisfy d-J^ = 0, = 0. Applying the general argument of the previous section, 
we conclude that we may choose the solution so that Eqs. |9.19| hold, so that JT" = jf = 0. 
This gives, altogether, 

^e^''d,^ = -ijri^, (10.25) 

so that the topological and Noether currents are proportional. As discussed in section ^ 
that is a consequence of the fact that E±2 satisfies ( p.l2| ). 

Let us turn to the Noether charge which here is simply the fermion number. It should 
be defined such that it satisfies the Poisson bracket relation 

i {'ip, QNoetherjp.B. = ^ (10.26) 

Since the coupling constant k was taken as an overall factor, this is satisfied by 

QNoether = k dx^J-f^il (10.27) 



39 



so that 

1 

Qtopol. = -^QNoether (10.28) 

As argued in general, this means that k should only take discrete values as expected, since 
our actions are related with the one of WZNW. 

Let us now construct the soliton solutions. The operators E±2 given in (|10.8| ), lie in 



the homogeneous Heisenberg subalgebra generated by if", with the commutation relations 
( |10.1| ). Such a subalgebra has no generators of grade ±1 for the principal gradation. There- 



fore, the model under consideration has no vacuum solutions of type ( |2.20|) . Then, from 
(13), we get 

S± = E±2 = mH^\ (10.29) 
We perform the dressing transformation starting from the vacuum solution ( p. 18 ), namely 



ip = r] = ip = ip = 0, ; z/ = -^mlx+x- = z/q. (10.30) 

Now, let I Ao ) and | Ai ) be the highest weight states of two fundamental representations of 
the affine Kac-Moody algebra sl{2)^^\ respectively the scalar and spinor ones. Then, from 
( p.l6| ) with 1] = 0, we obtain the solutions on the orbit of the vacuum ( |10.30| ), 



' ^ ' e-(^-^«) = (Ao I G I Ao), 

(Ao I G I Ao)' 



[m (\o \ ElG \ Xo) ~ fm (Xi \ ElG \ Xi) 

Vr = \ - TT — , ^ , ? , , i'R 



where 



(Ao I G I Ao) V 2 (Ai I G I Ai) 

/m (Ai I GE^_ I Ai) ~ /m (Ao | G E~^ \ Xp) 

IpL = —\- ■,WL = —\-. — , (10.31) 

V ^ (Ai I G I Ai) ' V ^ (Ao I G I Ao) 



G = e^+^+ e"^-^- pe^-^- 6"^+^+. (10.32) 



In order to get the soliton solutions, we choose the fixed mapping p to be an exponentiation 
of an eigenvector of £± (solitonic specialization); namely, p = e^, with [S± , V] = uj±V. 
Therefore, 

G = exp (^e^ with F = cx;+x+ — oo-X- = 'j {x — vt) . (10.33) 

In this case the eigenvectors of £± are 

V±{z) = ^ z-'^El. (10.34) 
ne2 



Indeed, 



[S+,V±{z)] = ±2mzV±{z) =uj^V±{z), (10.35) 
[g_,V±{z)] = ± V±{z)=u^V±{z). (10.36) 



z 
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The solution, associated with V+{z), is 



u = uo, = = 0, V 
while those, associated with V^{z), is given by 
Z/ = Z/Q, y9 = ^ = 0, ij = 

where 

1 



z 
-1 



(10.37) 



1 

1/z 



2m{zx+ x„) = 7 (x — ft) 



(10.38) 



;i0.39) 



The masses of these solutions are obtained from ( [7.17| ). Here the relevant state | Ag') in 
(Ol) is 



I K') =1 Ao)® I Ai). 
Using level one vertex operators, one can verify that 

{\i I (V±(2))" I A,) = , forn > 1 and i 



0,1. 



;io.4o) 



;io.4i) 



Therefore, Ny = in ( |7.1(j| ), and from ( [7. IT] ) one gets that the masses of the solutions 
( |10.37| ) and ( p.0.38| ) vanish. Such solutions correspond to the objects which travel with 
velocities v = ±(1 — 2;^)/ (1 + ^;^); and keeping z"^ > 0, one has | f |< 1. Therefore, 
these solutions cannot be interpreted as solitons (particles), since they would correspond to 
massless particles traveling with velocity smaller that light velocity. We should interpret 
them as vacuum configurations, since they have the same energy as vacuum ( |10.30| ). 

The true soliton solutions of the system are constructed as follows. Notice that V-^{z) 
and V^{—z) have the same eigenvalues. Therefore, any linear combination of them, leads to 
solutions traveling with a constant velocity without dispersion. So, we let 



V{a±,z) = Vtia+V+iz)+a-V-i-z)); 



[£+ , V{a±, z)] = 2mz V{a±, z) 



[S^,V{ai,z)] 



2m 



V{a±,z), 



(10.42) 
(10.43) 



and so cu^ = 2mz and tu. = The particular factor \/i is chosen such that the reality 
condition will be obeyed with a_ = a\. Again, using level one vertex operators, one can 
verify that|^ 

(As' I V{a±, zY I As') = for n > 4. (10.44) 

Therefore, Ny = 4 in (|7.16|) , and from ( |7.17| ) with ■^/'^ = 2, and z/' being the highest root of 
s/(2), one gets that the mass of such solutions is 

M = 8km = 2km^, (10.45) 

-'^''Notice that the truncation occurs for powers greater than 4, and not 2, because | As') hes in the tensor 
product representation, see ( 10.40| ) 
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where k is the couphng constant appearing in the Lagrangian ( 10.18| ), see ( |6.12| ). The 
solutions generated by (|10.42|) , have two parameters, namely a±. One parameter is always 
present, because one can scale an eigenvector of S± without changing the width 7 and velocity 
V of the soliton, obtained from the eigenvectors uj±] see ( |5.19|) . However, in this case, the 
second parameter comes from a symmetry. As we have pointed out in (|5.21|) , associated to 
the fixed element p = e 



V{a±,z) 



we have an orbit of equivalent solutions due to the global 



transformations ( |10.21|) , 



V{a±,z) ^ Vi [a+ e''' V+{z) + e-''^ V-{-z)) . (10.46) 
The explicit form of the solutions generated by ( |10.42| ), is obtained using (|10.31|) , 



V 



log 



2r^ 



tae 



2r 



log ( 1 + iae 



2r 



-m: 



X- 



0; 



and 



V l-jfTe2r 



r / i-io-e^r 

me ri 



(10.47) 

(10.48) 
(10.49) 



(10.50) 



where F is given in ( |10.39| ), and a = a^a-z/4. Keeping m and z real, we have the mass M 
of the soliton, from (|10.45|) , real and positive, and also the parameters 7 and v (|10.39|) are 
real. The reality condition is obeyed if a_ = a+, as anticipated. At this point, it is useful 
to re-express the expressions just given in terms of the physical parameters of the soliton. 
Using equations ( ^.231) , and ( |10.12[) , one deduces that 



7 = / Vl-v^ , z = \J{1- v)/{l + v). 
Moreover, since a± are complex conjugate, we may write 



The dependence upon space-time appears only through A/o-exp(r). We will write 

\fae^ = exp((7(a; — xq — vt)) 
where xq is the position of the soliton at time zero. Then we have 
(fi = 2i arctan (^exp (2m^ {x — Xq — vt) / \J\ — v' 
which is the sine-Gordon soliton. The Dirac fields are given by 

1 



l+v) 

1/A 



l+ze 



V 



(10.51) 

(10.52) 

(10.53) 
(10.54) 

(10.55) 



l—ie 



2m^{x — XQ—vt)/y/\ — v'^ ' J 



^by convention, we choose a to be positive 
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and if) is the complex conjugate of %p. Thus the only parameters are the soliton mass and 
velocity, together with the angle 9 which reflects the global invariance ( |10.21|) . Notice that 



the sign of the tolopogical charge can be reversed by reversing the sign of z. Therefore, the 
solutions ( p.0.54| )- (|10.55|) contain the sine-Gordon soliton and anti-soliton. 

Finally, we come to the very important feature of the present model already mentioned 
above in general, namely it is clear from the explicit expressions Eqs. |lQ.55 that ip vanishes 
exponentially when x — Xq ^ ±cxd, so that the Dirac field is confined inside the soliton. 
That this must be true is of course a general consequence of Eq. |10.25| which may be verified 
directly on the explicit solution. This phenomenon has been much studied for electron 
phonon systems. Models of a similar type describe the electron self-localization in quasi-one- 
dimensional dielectrics (for recent reviews see |32|, [Q). At low temperature these systems 
go over to dielectric states characterized by charge density waves which can be constructed 
on the basis of the Peierls model. The continuous limits are described by Lagrangians similar 
to Eq. |10.18| . Discussing this important issue is beyond the scope of the present article, so we 
will not dwell upon it here. Let us simply recall that the typical example of the polyacteline 
molecule was much discussed in connection with fermion number fractionization [0, 
Clearly, on the other hand one may regard our soliton solution a sort of one dimensional bag 
model for QCD. In this connection let us note that, if we introduce the two-by-two matrix 
exp(?7 + 2(^975), we may rewrite the Lagrangian Eq. |10.18 



U 



as 



— (tr 
16^ 



M 2 



1 

- -tr 
2 



U 



tr 



u- 



m. 



;i0.56) 



which is similar to a two-dimensional version of the low energy effective action for QCD (see 
e.g. H). 



10.2 Case / = 3 

In this case we choose 

and so 



E_3 = g+^;' + g-^l'; 

[E^,E_^]=?>q+q_C = f3C. 



Introduce also the notations 
El 



q+El + , = q+E-^ + 

h = q+El-q^E'_, f^i = q^E-'-q^E'_ 



(10.57) 
(10.58) 



;i0.59) 



The fields of the model are those introduced in (|10.7|) , and the matter fields 'i/'^/^, i = 1,2, 
and x± defined as 



^1+ = ^Ifi + X+Ei, 



= i^Rf2; 



(10.60) 
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According to ( ^.201 ), V'Ij/L' = 
denote by ip^. 

One can easily verify that 



1,2, are the components of two Dirac fields which we shall 



[E3 , [-E_3 , E^i] 
[E3 , [E.3 , H'] 

[E3, [E^3,fi] 







±1,±2. 



Therefore, from ( [7 .51 ), the masses of the particles are 



nir. 



m 







m^r = 4:^q+q- , z = 1,2 



and so we have to choose q±, such that q+q- > 0. 

Then, from ( ^■13| )-( p.l(j| ), the equations of motion are 



d+d- (f 



q+Q- (( 
+ 

(( 



2(p 



2X+ - i^RX 
2 2 



2 V 

1 



R 



(< 
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-2</3 



X- (e- 



2y^ (^-e^V^j + 2^'^! 

d+^Pl = 
d-X+ = 
d+X- = 

Notice that these equations are invariant under the CP transformation 



-2^ 



./q^e^^lj'n (x- (e^^ - e'^^) - (. 
0. 



1 e^'P + 



;io.6i) 



(10.62) 



X_ , ^ > ^li ^ -^L . X+ ^ X- , <^ ^ <^ ^ 



(10.63) 

(10.64) 
(10.65) 
(10.66) 

(10.67) 

(10.68) 

(10.69) 

(10.70) 
(10.71) 



7] ^7]. (10.72) 



Since, there are no generators of Qq, given in (|10.6| ), that commute with E^^^ in (|10.57|) , we 
do not have any gauge symmetry of the type ( p.29| )-( p.31| ). In the linear approximation for 
the V' and x fields with 93 = and r] = constant, one verifies that the ip field equations may 
be deduced from a Lagrangian of the type Eq. ^.27| , with ip^ ip, and ip'^ ^ ip. Such is not 
the case for the x fields although they are massless in the linear approximation. Indeed, x± 
has weights (|,0) and (0, |) respectively, so that we cannot introduce a conjugate field to 
write down a covariant kinetic term. As a matter of fact, we have been unable to derive the 
above field equations from a local action. 
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The sine-Gordon (or sinh-Gordon) model is a submodel of the system ( |10.63| )-( |10.71| ). 
Indeed, = x± = ''7 = is a solution of the equations of motion, and then iip {ip) has to 
satisfy the sine-Gordon (sinh-Gordon) equation. 

The operators i^-ta belong to the principal Heisenberg subalgebra of sl{2Y^^ 

[E2m+1 , E2n+l] = q+q-{2m + l)C 5n+m+lfi (10.73) 

where 

^2m+i = "^eZ (10.74) 

They are eigenvectors of the grading operator ( |10.5| ) 

[Q , E2™+l] = (2m + 1) E2rn+l (10.75) 

The adjoint action of E2m+i is diagonalized by the operators 

OO OO -| 

V{z) ^ -Vq^ z-^-H-+ z-'--' {q+El-q.E!l-^')--,/q^C 

n=— OO n=— OO 

OO 

= E ■2""^n. (10.76) 



where 

Indeed, one gets 



[Q,Vn]=nVn (10.77) 



[E2m+i , Viz)] = 2^q^ z^^+^ V{z) (10.78) 
Therefore, V{z) are eigenvectors of E±-i with eigenvalues 

uj± = 2^fq^ z^'^ (10.79) 

Notice that, if one shifts z — >■ loz^ with lJ^ = 1, the eigenvalues do not change. Therefore 

V{aj, z) = aQV{z) + aiV{uz) + a2V{u'^z) (10.80) 

are also eigenvectors of -£^±3. Therefore, if one performs the dressing transformations from 
the vacuum (|2.18|) , namely 

(/9 = ?7 = X± = '^"^ = V'^ = 0; = ~3g+g_ z+x„ (10.81) 

one obtains solutions traveling with constant velocity, without dispersion, by taking the 
constant group element p in ( ^.16|) as the exponentiation of (|10.80|) (see (|5.19| )). 

The operator V{z) introduced in ( |10.7(j| ) can be realized through the principal vertex 
operator construction |l^]. Using such contruction, and taking into account that the 
highest weight state ( |10.4U| ) lies in the tensor product representation, one gets 

for n > 2, with just one non vanishing a/s (10.82) 
for n > 4, with two non vanishing a^'s (10.83) 
for n > 6, with all three a/s, non vanishin^lO.84) 



(As' 




z)T 


1 As') 


= 0; 


(As' 


(F(a„ 


z)r 


1 As') 


= 0; 


(As' 


(V(a,, 


z)r 


1 As') 


= 0; 
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Using the mass formula ( |7.17] ), one gets that the masses of the sohtons created by the 
operator V{aj,z) in ( p.0.8CI| ), from the vacuum ( |10.81| ) are 



Ml = -k ^q^q^ = 2k ; with just one non vanishing a/s (10.85) 
16 

M2 = — k y/q+q- = Akra^p ; with two non vanishing Oj's (10.86) 
o 

M3 = 8 A; y/q+q^ = 6km^ ; with all three a/s non vanishing (10.87) 
Now, if one performs the dressing transformations from the vacumm (|2:20| ), namely 



ip = r] = ijj^ = ijj'^ = 0, X± = ^ v = —"iqj^q- (10.88) 
with = const., then the soliton type solutions are created by eigenvectors of 

£± = E±3 + c^^ii (10.89) 
Again, V{z) are eigenvectors of £± with eigenvalues 



Therefore from ( [7.17| ) one gets that the mass of the corresponding soliton is {Ny = 2 from 

mm) 



M' = -k ^q+q_ (1 + ctcT + c{z^ + cl / z^) (10.91) 
o 

The spectrum of the soliton solutions of such model is rather complicated. There re- 
mains to perform a more detailed analysis of those solutions, in order to understand, among 
other things, the physical consequences of the existence of (at least) two classical vacuum 
configurations. 

11 Outlook 

The models we introduced in this paper constitute a generalization of the affine Toda systems, 
in the sense that they contain matter fields coupled to the usual (gauge) Toda fields. We 
believe such models open the study of a variety of massive integrable theories with very 
interesting physical properties. As we have already mentioned, from the point of view of 
non perturbative aspects of quantum field theories, we hope these models will be useful, 
as a laboratory, in the understanding of the quantum theory of solitons, some confinement 
mechanisms and also to obtain exact results on the strong coupling regime. On the other 
hand, these systems can be used to describe very interesting phenomena in condensed matter 
physics and statistical mechanics like electron-phonon systems, electron self-localization in 
quasi-one-dimensional dieletrics and polyacetylene molecules. 

The next step in achieving such program is to consider the quantum theory of these 
models, trying to obtain exact results by exploring some of their special physical proper- 
ties. In this sense the most promising class of models is that described in section P, which 



simplest example corresponds to the model associated to s/(2), described in subsection |10.1 
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That class of models possesses a U{1) Noether current depending only on the matter fields 
which, under a special gauge fixing of the conformal symmetry, is equivalent to a topologi- 
cal current depending only on the (gauge) Toda fields. We believe the full consequences of 
such equivalence have not been understood yet and, the use of it, may shed light on several 
non perturbative aspects of the theory. We have already pointed out that such equivalence 
leads, at the classical level, to the localization of the matter fields inside the soliton. It is 
very plausible that, at the quantum level, this gives a confinement mechanism, which can be 
regarded as a one dimensional bag model for QCD. Another special aspect of these models 
concerns the mass formula. The masses of solitons and particles are both determined by the 
eigenvalues of the operators E±i appearing in the fiat connection defining the models. In 
addition, the soliton masses have a topological character due to the spontaneous breakdown 
of the conformal symmetry. The models presenting the equivalence between Noether and 
topological currents have an additional feature; the masses of particles and solitons are pro- 
portional to the f/(l) Noether charge. This situation resembles very much the one of four 
dimensional gauge theories with Higgs in the adjoint representation and in the BPS limit. 
These facts may indicate the existence of a sort of duality in these models involving solitons 
and particles [0. These points certainly deserve further study and we hope to explore 
them in a future work. They may have a direct connection with super symmetric Yang- Mills 
theories along the line of ref.|I| 

In addition to those, there are still many interesting practical problems to be solved 
following the lines of the present paper and references herein. We mention first the study 
the VT-symmetries of the system under consideration, in particular, generalizations of the 
VT-algebras of the standard Toda systems, ly-geometries associated with them in the spirit 
of [^, [^], [Q, including differential and algebraic geometry setting of the Toda systems 



coupled to the matter fields. For this questions, the formulation of the corresponding problem 
in terms of Lax operators of a Generalized mkdV hierarchy [^| should be useful. Second, 
it is very believable that the general solution for the matter fields can be presented in a 
compact form; in other words, it seems to be possible to resolve the recurrent formula ( 3.13 ). 
Third, it would be useful to obtain, in some suitable parametrisation for the mappings B and 
F^, an explicit expression for the Lagrange function or the effective action corresponding to 
our system (|2.7| )-( |2!9| ). Having such an expression, one can directly get the formula for the 
energy-momentum tensor, and hence calculate in a more simple way the masses of solitons 
generated by the spontaneously breakdown of the conformal symmetry. In cases where this 
is not possible, one can work with the WZNW fields, in terms of which the matter fields can 
be written locally. 

The construction of the multisoliton solutions can be made in a straightforward way 
using the methods presented in this paper. That would be important in the study of the 
classical scaterring of the solitons, which could give us valuable information to construct the 
corresponding S-matrix, using a bootstrap programme |^T], Here, the question of 

massless particles is of special importance |^ . 



Another point to be explored is the construction of the local conserved charges of these 



systems using for instance, the methods of refs. P5| , where the fiat connection is gauge 
transformed into an abelian (Heisenberg) subalgebra of the affine Kac-Moody algebra Q. 
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